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PREFACE 

The purpose of this manual is to provide a rapid reference 
for the CALCOMP Plotter users at TSC. The present system uses 
FORTRAN 10 on the DECsystem-10 to generate internal files to 
be plotted in STAND-ALONE Mode on the CALCOMP Model 763. 

An extensive software capability is being planned and has, 
in part, been implemented for use in a wide variety of current 
data processing applications. Routines are available that 
construct and annotate axes and grids with symbol types and 
alphameric lettering of many sizes and angles of printing. 
Numerous routines are currently operational for two dimensional 
display of data sets (both scaled and edited) . In the future, 
we plan to add routines which will provide a 3 dimensional 
plotter capability. Curves may be constructed through the data 
and specified points identified with special symbols 
automatically. Information may be filtered and smoothed to 
reject extraneous noise if required. Programs are also avail- 
able to interpolate between available points and establish 
missing points, or to extrapolate from various inputs to extend 
the curve. These routines are found in the SSP Library of 
programs and are FORTRAN callable. 

The full utility of the system depends upon the use by the 
applications programmer. Since the entire system is based on 
simple FORTRAN type call statements, a user may incorporate 
the plotting capability directly into a data reduction or 
analysis program. It is the users choice as to whether to 
generate his plot using the basic standard CALCOMP software or 
to generate the plot using the individual routines. 
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Some of the information in this manual is adapted from the 
California Computer Products, Inc., literature. Other segments 
are from the Plotting Manual which was generated by the Analysis 
and Computation Division at the NASA Langley Research Center. 
Still other portions of this manual are written by Messers. T. 
Talbot and R. Shah of KHL under Support Services Contract. 
These two gentlemen were instrvmiental in modifying (for the 
DECsystem-10) the majority of the software presented in this 
doctiment. Appropriate acknowledgements as to the source of the 
particular routine will be found with the routine write-up. 
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INTRODUCTION 



The pxirpose of this manual is to provide a common 
reference for those ADP users at the Transportation Systems 
Center who are engaged in the development and use of Computer 
Plotting software. Included are general descriptions of the 
currently cataloged routines; procedural guidelines for pre- 
paring, documenting and processing the plotting routines; and 
a reference manual of program write-ups with examples. 

This manual is prepared and maintained by the Computer 
Services Division. A copy will be issued to a user by a 
request to the PIC (Program Information Center) , Mr, M, 
Compagna, PIC Coordinator. Updates are mailed to each manual 
holder as the situation demands, and it is the holder's 
responsibility to incorporate them into his manual binder. 

Information of a temporary or announcement nature will 
be mailed out as needed to the user. 

Questions on the contents of this manual should be 
addressed to the ADP Applications Branch (862) . Questions 
on distribution should be directed to the PIC (862). 



ROUTINE IDENTIFICATION 



A plotting routine, identified by a Program Number , will be 
assigned when the program is accepted as part of the CALCOMP 
Library of programs . A program number has the form of one letter 
followed by five digits. The letter, "S" in this case, desig- 
nates the basic category of the program, i.e., PLOTTING; the next 
two digits, 30 indicate that this is CALCOMP plotting vs. 31 
which is SC4020 Plotting and 39 which is other Plotting; the 
remaining three digits are assigned in the following manner by 
subject category: 

000 - 099 Line applications; line types, labeling, etc. 

100 - 199 Scale applications 

200 - 299 Axis applications; grid types, etc. 

300 - 399 Special symbols; alphabets, etc. 

400 - 499 Special functions; curve smoothing* french c.urv^, 
etc. 

500 - 599 Special forms; vibrating spring, etc. 

600 - 699 Pictoral applications; maps, shading, etc. 

700 - 799 Unique plots; 3D applications, etc. 

800 - 899 Specialty routines; card deck, etc. 

900 - 999 Miscellaneous 

The plotting routines are identified by name and are cata- 
loged according to their function. Either the name or the number 
will find the routine. 
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ROUTINE DOCUMENTATION 

Routines for which documentation will be prepared are 
expected to be used as a closed loop within other programs and 
will have general interest and application. The following out- 
line will be followed in preparing routine documentation. 

Title: - The title includes the words PROGRAM, SUB- 
ROUTINE or FUNCTION and the code name; e.g., SUBTROUTINE 
EGNC . 

Language; - Indicate FORTRAN, or other appropriate language. 

Purpose ; - A brief statement of what the routine does. 

Use: - Either control card for a program, function state- 
ment for a function, or calling sequence for a subroutine. 
All variables will be clearly defined. Diagnostic or error 
messages will be listed. 

Restrictions ; - Type, dimensions and range of arguments, 
and special hardware requirements. 

Method ; - Derivation, equations, or other explanation of 
the method used. 

Accuracy : - A statement of accuracy or "not available" 
statement. 

References ; - Published material related to the material 
used. 

Storage ; - Storage requirements (words) . 

Subprograms Used; - List of standard subroutines with 
number of words required. 
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other Coding Information ; - This should be included when 
pertinent and may include such information as timing, 
explanation of error messages or diagnostic statements, 
special subprograms, etc. 

Source ; - Organization and person's name (if not developed 
at TSC) . 

Responsible Person; - Name of programmer responsible for 
providing the program for inclusion into the library. 
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ROUTINE UTILIZATION PROCEDURES 

The following are descriptions of the different available 
plotter routines. The program number and title of each routine 
are followed by necessary information for its use. Note that 
in the argument list of each routines, FORTRAN convention are 
maintained. Any argument starting with (I, J, K, L, M, N) is a 
fixed-point integer value; all others are floating point numbers 
with the exception of ALPHA arrays. 

Disclaimer : Although each routine has been tested by its 
contributor, Data Services gives no warranty as to the proper 
functioning or accuracy of the routine. No responsibility is 
assumed in the distribution of these routines by Data Services. 

A list of the presently available routines is tabulated 
on the following pages. This list includes the CALCOMP basis 
software which is compressed of six routines that can be called 
by your FORTRAN program. They, in turn, call upon each other for 
certain operations. A particular application may not need all of 
these subroutines, depending on the user's requirements. The six 
basic routines are as follows: 

PLOT* converts all pen movement specifications from 
inches to actual plotter commands, and outputs 
these to the attached device. 

SYMBOL* draws any sequence of alphanumeric characters, 
as well as special point-plotting symbols. It 
calls only PLOT. 
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NUMBER* draws the fixed decimal equivalent of an 
internal floating point nxamber. It calls 
only SYMBOL, which then calls PLOT. 

SCALE* examines a data array to determine an 
optimtun starting value and a scaling 
factor for use by AXIS and LINE in con- 
verting data units to plotter page dimensions. 
It is the only subroutine that does not call 
any other because it does no actual plotting. 

AXIS* draws an axis line with the appropriate scale 
annotation and title. It calls SYMBOL and 
NUMBER as well as PLOT. 

LINE* plots a series of scaled data points defined 

by two arrays (X and Y) , connecting the points 
with straight lines, if desired. It may call 
SYMBOL as well as PLOT. 
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ROUTINE LIST WITH PROGRAM NUMBER 
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LINE APPLICATIONS 



S30000 
S30001 
S30002 
S30003 
S30004 
S30005 
S30006 
S30007 
S30008 
S30009 
S30010 



PLOTS* 

PLOT* 

PLOTA 

LINE* 

ALINE 

CLINE 

NOMEN 



♦CALCOMP BASIC SOFTWARE 
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SCALE APPLICATIONS 



S30100 
S30101 
S30102 
S30103 
S30104 
S30105 
S30106 
S30107 
S30108 
S30109 
S30110 



FACTOR* 

SCALE* 

ASCALE 

LSCALE 

LOGSCA 



*CALCOMP BASIC SOFTWARE 
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AXIS APPLICATIONS 



S30200 
S30201 
S30202 
S30203 
S30204 
S30205 
S30206 
S30207 
S30208 
S30209 
S30210 
S30211 
S30212 
S30212 



AXIS* 

AXIS2 

AXIS3 

LAXIS 

AXLABL 

LOGAX 

DIMENS 

PAXIS 



GRID 

BORDER 

LIGRID 



*CALCOMP BASIC SOFTWARE 
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SPECIAL SYMBOLS 



530300 SYMBOL* 

530301 NUMBER* 

530302 ENUMBER 

530303 LABEL 

530304 NOTATE 

530305 PARROW 
S30306 

S30307 
S30308 
S30309 
S30310 



*CALCOMP BASIC SOFTWARE 
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SPECIAL FUNCTIONS 

S30400 
S30401 
S304Q2 
S30403 
S30404 
S30405 
S30406 
S30407 
S30408 
S30409 
S30410 
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SPECIAL FORMS 



S30500 


ARC 


S30501 


BLIPS 


S30502 




S30503 




S30504 




S30505 




S30506 




S30507 




S30508 




S30509 




S30510 
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PICTORAL APPLICATIONS 



S30600 


WORLD 


S30601 


SHADE 


S30602 


SHADIT 


S30603 




S30604 




S30605 




S30606 




S30607 




S30608 




S30609 




S30610 
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UNIQUE PLOTS 



S30700 THREE 

S30701 

S30702 

S30703 

S30704 

S30705 

S30706 

S30707 

S30708 

S30709 

S30710 
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SPECIALTY ROUTINES 



S30800 
S30801 
S30802 
S30803 
S30804 
S30805 
S30806 
S30807 
S30808 
S30809 
S30810 



DECKS 
DOCPRT 
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MISCELLANEOUS 



S30900 



f\ T '^ 



ROUTINE DESCRIPTIONS 
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LINE APPLICATIONS 
S300XX 
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IDENTIFICATION : PLOTS ' 



Purpose : 



Restrictions: 



S30000 

The call to this subroutine initializes 
the Plot subroutine. The buffer areas 
are established and all devices speci- 
fied. 

This routine is called only once before 
a call to any other plotting routine is 
made . 



Usage: 

Calling Sequence: 

Storage Required! 



Method : 



Call PLOTS (0,0,16) 

16 - Logical unit number for plotting 
output. 

This routine is one of the standard 
CALCOMP routines and can be found in 
the PROGRAMMING CALCOMP PEN PLOTTERS 
manual from CALCOMP, Inc. 



Original Developer 



.* 



CALCOMP, Inc. 
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IDENTIFICATION: 



PLOT* S30001 



Purpose : 



This routine moves the pen either in an 
up or down state to a new position. 



Restrictions : 



None 



Usage: 
Calling Sequence; 



CALL PLOT CX, Y, ± X) 

Y, X - The terminal position coordinates to 
which the pen is moved, in inches, 
from the current reference point. 

32 ■ Controls pen up/ down status, and 
origin definition. 

JC. ■ 2 - pen is down during movement 
drawing a line. 

il = 3 - Pen is up during movement. 

■J^=-2-3-A new origin is defined. 



Storage Required: 



Method: 



Original Developer: 



This routine is one of the standard 
CALCOMP routines and can be found in the 
PROGRAM^UNG CALCOMP PEN PLOTTERS manual 
from CALCOMP, Inc. 
CALCOP, Inc. 



lie 1 



Where, 



Identification; Subroutine PLOTA S30002 

Purpose; 

To draw a figure or curve rotated about a point (XW, YW) 

Usage; FORTRAN Calling Sequence 

CALL PLOTA (X Y, I) 

Additional input is through common block lANG 

COMMON/IANG/ SINAN , COSAN , XW , YW 



X, Y and I are exactly the same as in the CALCOMP sub- 
routine PLOT. 



SINAN 
COSAN 
XW, YW 



Sine of angle of rotation 
Cosine of angle of rotation 
Point of rotation 



Other Subroutine Required ; 
Storage Required: 
Original Source 
In-house Submitter: 
In-house Contact: 



PLOT 

DECsystem (61) g Locations 

Tom Talbot/KHL 

Tom Talbot/KHL 

Tom Talbot/KHL 
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r 

r A TFST PROGRAM TO <;HOW THE USF OF SUBROUTINES PLOTA AND ARC 

r 

COMMON/IANfi/«^INAN,ro<;AN«XW,YW 

DIMFN5.I0N XCJ) .Yf?) ,TPFN(5) 
r 

r «;ft data to draw an arrow 
r 

DATA X/n,,.«5,.A,.5,,4/ 
OATA v/n,»n^,^l ♦^.♦-.1/ 
DATA TPFM/?,?,?,-^,?/ 

r 

CALL PLOTS (n,n, 16) 
r 
r <--FT ORTC^TN 

r 

CALL PLOT(0,»-n.,-3) 
CALL PLOTC^.^A.,-^) 

r 

r SFT POINT OF ROTATION 

r 

xw=-?.o 

r 

r «:PT AMf;| r 

C 

no 1 1=1,1? 

ANG=FLOATf I1»3(^.»0,0174S3'^ 

5;iNAN=SIN(AMG) 

COSAN=CO<;(ANG) 

r 

r DRAW ARROWS (PLOTA -ROUT INF) 

r 

no 1 j=i,5 

CALL PLOTA(X{J),Y{J),TPFN(J)) 
1 rONTTMUF 

r 

r DRAW A riRCLF 
r 

CALL ARr(XW,YW,n,,'?60.,2,5,3) 
C 
r END PLOT 

r 

CALL PLOT(8,fO,,-3) 
STOP 

FMD 
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IDENTIFICATION: 



LINE 



S30003 



Purpose : 



Restrictions: 



Usage : 
Calling Sequence: 



This routine produces a line plot of 
the pairs of data values in two arrays . 
The dot points may be represented by 
centered symbols and/or connecting 
lines between points. 

Scaling parameters corresponding to FV 
and DV must immediately follow each 
array. If these parameters have not 
been computed by the SCALE subroutine, 
they must be supplied by the user. 

Call LINE (X, Y, N, INK, + LINETY, ISYM) 

X - Array of abscess a values and 
scaling parameters C+2) for the X-array. 

Y - Array of ordinate values and scaling 
parameters (+2) for the Y-array. 

N - Number of points to be plotted in 
array X. Y-array must have same number 
of points. 

INK - Increment to use in gathering data 
from X and Y arrays. Must be same as 
from SCALE routine or be a calculated 
value given by the programmer. 



2.15.3 



+_ LINETY ' Describes the line type to 
be drawn through the data points. The 
magnitude determines the frequency of 
plotted symbols. 

NOTE: = = > No symbols plotted. 
= +K = > points straight line 

connected. 

= -K = > only symbols plotted; no 

connecting lines. 

ISYM - Value from through 13 for 
symbol type to be used. 



Storage Required : 
Method: 



This routine is one of the standard 
CALCCM* routines and can be found in the 
PROGRANMING jCALCGMP PEN PLOTTERS manual 
from CALCCWP, Inc. 



Original Developer; 



CALCCMP, Inc. 



2.15.3.1 



v^here , 



Identification; Subroutine ALINE S30004 

Pur pose ; 

To drav; line through points (Xi, Yj^l 

Usage; FORTRAN Calling Sequence 

CALL ALINE Dw Y/ N, XM, DELX, YH, DELY) 

X, Y Array of coordinates to be plotted 

N Number of points 

XM, YM Minimum X and Y respectively on the axis 

DELX, Unitr; per inch on the X and Y axis 
DELY respectively 

ComiTientc : 



This routine draws a line betv/een points,. and the 
scale factors need not be included in the arrays as is 
required by Calcomp routine LINE, 
other Subroutines Required: PLOT 

Storage Required: DECsystem-lO U.02)g Locations 

Original Source: Tom Talbot/KHL 

In-house Submitter: Tom Talbot/KIIL 

In-house contact: Tora Talbot/KHL 



2.15.4 ^ 



Identification; Subroutine CLINE s 3000 5 
Purpose ; 

To draw different types of lines connecting 
points (X. , Y. ) . 

Usage; FORTRAN Calling Sequence 
CALL CLINE (X, Y, N, XM, DELX, YM, DELY , NC) 
where, 

X,Y Array of coordinates to be plotted 

N Number of points to be plotted 

XM, YM Minimum X and Y respectively on the axis 

DELX, 

DELY Units per inch on the X and Y axis 
respectively on the axis 

WC Line type to be used. NC ranges from 
1 to 5. 

COMMENT ; 

CLINE connects points in the arrays X and Y 

in one of the five Line Types, as follows; 



LINE FOR NC=1 

. -LINE FOft.NC^a 

LINE FOR NC33 

LINE FOR NCrq 

-^-LINE FOR NCrS 

2.15.*5 



other Subroutines Required; 
Storage Required : 

« 

Restrictions: 
Original Source: 
In--house Submitter: 
In-liouse Contact: 



PLOT 

DECsystem-lO (311) g Locations 

NC ranges from 1 to 5 

Tom Talbot/KIIL 

Tom Talbot/KJIL 

Tom Talbot/KJIL 



2.15.5.1 



Identification; Siibroutine NC»IEN S30006 

Purpose ; 

To draw an inch of line type (as in CLINE) followed 
by a title for annotation. 

Usage ; FORTRAN Calling Sequence 
CALL NOMEN (X, Y, NC, HT, BCD, NCHAR) 



where , 



X, Y Start point of title in inches 
NC Line type (1-5) , to be used with title 
HT Height of characters to be used in title 
BCD Array containing characters for title 
(H format accepted) 

NCHAR NxMtiber of characters in title 



t(\ in 



cU, 



Other Subroutines Required: 
Storage Required: 
Original Source: 
In-house Submitter: 
In-house Contact: 



PLOT, CLINE 

DECsystem-10 (110) g Locations 

Tom Talbot/KHL 

Tom Talbot/KHL 

Tom Talbot/KHL 
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r 

r A TF9T PROr^RAM TO .<;HOW THF USF OF <;i IBROUTINES ALINE ♦AXIS2 ,CLI NF , 
r AND NOMFN 



r 



DIMFN,<;ION XfinO) ,Y(100) ♦YC(lOn) ,YE( ion) 

r 

CALL PLOTS(n,0»16) 

r 

r «;ft opTfiTN 
r 

CALL PLOT(n,,-n ,,-'^) 
CALL PLOT(n,'5,0,«5,-3) 

r 

r genfratf data for cijrvf*; 
c 

do 1 1=1,100 

X( T)eFL0AT(I)*.n5 

Y(n=.<;iN(x(n ) 
YC(i)=co^(X( in 
YE(n=Exp(Y(n ) 

CONTTNIJF 

r 

C EACH CALL TO AXIS? RETURNS THE APPROPRIATE SCALING FACTOR 

r X<;C AND YSC 

r 

CALL AXTS?(n. ,o,,5.»n.,l,,-6.,8HX - AXI«: ♦8,-l ^XSC) 

CALL AXIS2(0,,0,,2,,-2.,1,,4,,8HY - AXIS » 8,-l , YSC) 
C 

C DRAW A SINE CURVE WITH SOLID LINE —USE OF SUBROUTINE ALINE 
C 

CALL ALINE(X»Y.100,0,,XSC,-2.,YSC) 
C 

r DRAW A COSINE CURVE WITH DASH LINE 

C DRAW AN EXPONENTIAL CURVE WITH LONG-SHORT DASH LINE 
r USE OF StlRRDUTINF CLTNF 
C 

CALL CLlNF(X»YC,lon,o, XSC,-?. ♦ YSC, ? ) 

CALL CLINE(X»YF,100,0.,XSC,-2.,YSC,3) 
r 

C SEPARATE ANNOTATION ON CURVES — USE OF SUBROUTINE NOMEN 

r 

CALL NOMEN(3.»4,25,l,n,12»10HY = SIN(X),10) 

CALL NOMEN(3«,A,0,2,0.1?,10HY = COS(X).IO) 

CALL NOMFN(3.,3.75,3,n.l2,15HY » FXP( SIN ( X ) ) .15) 

r 

r FND PLOT 

r 

CALL PLOT(R,,0,,-'^) 

STOP 
END 



SIN(X) 
COS(X) 
EXP(SIN(X)) 




- AXIS 
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SCALE APPLICATIONS 
S301XX 
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IDENTIFICATION: 



FACTOR 



* S 



30100 



Purpose: 



This routine enables the user to enlarge 
or reduce the size of the entire plot 
by changing the number of plotter steps 
per inch of page coordinates. 



Restrictions: 

Usage: 
Calling Sequence: 



Storage Required: 
Method: 



None 

Call FACTOR (FAC) 

FAC - The ratio of new plot size to nor- 
mal plot size, e.g., if FAC = 2.0, the 
plot will be twice the size of a normal 
plot. 

This routine is one of the standard 
CALCCMP routines and can be found in the 
PROGRANMING CALCCMP PEN PLOTTERS manual 
from CALCCMP, Inc. 



Original Developer: * 



CALCCW, Inc. 



2.16.0 



IDENTIFICATIC»^: SCALE 



S30101 



Purpose : 



This routine is used to examine the data 
values in an array and to determine a 
starting value, either maximum or mini- 
mum, and a scaling factor. 



Restrictions: 



None: 



Usage : 

Calling Sequence: 



Note: 



Call SCALE ( A, AX, N, + INK) 

A - First location of array of points 

to be sealed. 

AX - Length of axis to which data will 
be scaled, must be greater than 1.0 
inch. 

N - Number of points in A-array. A should 
be dimensioned at least ACAR +2) 

+_ INK - Increment to select values. 
Usually 1. 

If positive, the A (N + 1) is a minimum 
and ACN +2) is positive. 

If negative, then A(N +1) is a maximum 
and A(N +2) is negative. 



2.16.1 



A(N + 1) is the starting value and 
ACN + 2) is the scaling factor. 



Storage Required: 



Method: 



This routine is one of the standard 
CALCCM' routines and can be found in 
the PROGPANMING CALCCMP PEN PLOTTER 
manual from CALCCM', Inc. 



Original Developer: 



GALC(M>, Inc. 



2.16.1.1 



scale : factor and stores it in X ( (N*K) +1+K) . The scale 
factor will be a power of 10 x; 2, 4, 5, or 10. The data 
in the array may be scaled to floating point inches by using 
a formula similar to the following: SV = (AE-MV)/SF, where 
SV = scaled value, AE = present value of array element, 
MV = either the minimum value or the value desired at the 
origin, and SF = the scale factor computed by the subroutine, 



Other Coding Information! 



Other Subroutines Required) 
Storage Required: 
Restrictions: 



Original Source: 
In-house Submitter: 
In-house Contact: 



Example: DIMENSION ORD (102), 
ABS (204) CALL ASCALE (ORD, 10., 
100,1,10.) CALL ASCALE (ABS, 15., 
100,2,10) 

None 

DECsystem-10 (401) g Locations 
The array must be dimensioned 
to include storage space for 
two extra elements per inter- 
leave factor. For example: 
N = 100, K = 1, DIMENSION X 
(102); N = 75, K = 3, DIMENSION 
X (231) 

Langley Research Center 
Tom Talbot/KHL 
Tom Talbot/KHL 



2.16.2.1 



Identification : Subroutine ASCALE S30102 
Purpose ; 

To compute a scaling factor for an array of numbers 
to be plotted and find the minimtmi data value within the 
array. 

Usage ; FORTRAN Calling Sequence 
CALL ASCALE (X, S, N, K, DV) 
where, 

X is the name of the array containing the 
floating piint values to be scaled. 

S is the length (floating point inches) over 
which the data are to be plotted (usually 
the length of one of the axes) . 

N is the number of data values in X from 

which points are to be plotted in accordance 
with K. 

K is the interleave factor which specifies the 
the sequence in which data are stored. 

= 1 indicates that values are stored sequentially 
= 2 indicates that values are stored in every 
other location in the array, etc, 
DV is the number of divisions per inch of the 

plotting paper to be used (should be: 10.0, 
20.0, 25,0, or 25.4), DV scales the polar 
magnitude array. 

Method ; 

This routine scans the elements in the array to find 
the minimum and maximum. ASCALE computes an adjusted minimum 
(origin value) and stores it in X ( (N*K) +1) and computes a 
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Identification : Subroutine LSCALE S30103 

Purpose ! 

To scale the array to be plotted with a log axis 
allowing user to specify the scaling of his data. 

Usage ; FORTRAN Calling Sequence 

CALL LSCALE (X, N, KVAL, ALONG, NCYCLE) 



where , 



X Array of values to be scaled 
N Number of points in the array 
KVAL Value of exponent (power of ten) at the first 
tic mark 

ALONG Length in inches of the log axis to be used 
NCYCLE Number of log cycles to be used 

Comments : 



LSCALE was written to allow the user to specify the 
scaling of his data. Unlike subroutine LOGSCA, the variables 
KVAL and NCYCLE are input by the user. The array X is altered 
as in subroutine LOGSCA. 



Other Subroutines Required: 
Storage Required: 
Original Source: 
In-house Submitter: 
In-house Contact: 



None 

DEC-10 

Tom Talbot/KHL 

Tom Talbot/KHL 

Tom Talbot/KHL 



(62) g Locations 



2.16.3 



Identifi cation; Subroutine LOGSCA 830104 

■ Ml ■ iii^ii.^ — I.I — I, — II ■ — ■■■■iM I ^11 IP m .M I . . m I III ■ i-i .III. »i — . ■ i» . m ill I. ■•- . 

Purpose ; 

To scale the array to be plotted with, a log axis and 
to return the variables required by the subroutine LAXIS. 

Usage; FORTRAN Calling Sequence 

CALL LOGSCA (ARRAY., N, ALONG, KVAL, NCYCLE) 

v;here, 

ARRAY Array to be scaled 

N Number of points in the array 

ALONG Length in inches of the log axis to be 
used 

KVAL Values returned bv LOGSCA for use by sub- 
NCYCLE r»dd ae^^«se^ 

Comment ; 

I This routine scales data to be plotted with a log 

axis and returns variables required by LAXIS. The values 

in the array are altered, so care must be taken by the user 

in its use after, the call. After scaling of both x and y 

arrays, the line may be drawn by subroutines ALINE or 

CLINE which do not require that the scaling factors be in 

ARRAY (N + 1) and ARRAY (N + 2) as does the Calcomp routine 

LINE. LINE may be used by setting ARRAY (N + 1) = 0. and 

ARRAY (N + 2) =1. The scale factors for M^INE or CLINE 

are: . 

XM=0. 

DELX =1. i 

■ . • i 
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other Subroutines Required: None 

Storage Required: DECsystem-10 (155) g Locations 

Original Source : Tom Talbot/KIiL 

In-house Submitter: Tom Talbot/IOIL 

In-house contact: Tom Talbot/KHL 
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c 

C A TEST PROGRAM TO SHOW THE USF OF SUBROUTINES ASCALE»LOGSCA tLAX IS 

r 

DTMFNSTON XAROAY fl V»FX ( 7) 

DATA FX/'>9fi.,2<»9.,no,,5<9.,2<»,»0.»0,/ 

CALL PLOTSf'^fOtlfi) 
r 

r e^FT ooiriN 
r 

CALL PLnT(n.,-ll.,-3) 

CALL PLOTH. ♦!.♦-'?) 

r 

r SFT srALF FACTOR 

c 

CALL ASCALE|FX»7,,5,1«10,) 

r 

r PLOT AXT<; 

r 

CALL AXIS?('^.,n,,,l,n.,,02.-«5,»-97HPROBARILITY OF SUCCESS AFTFR N 
»TRIALS,37,^,DFLX) 

r 

r Fwn PLOT 

r 

CALL PLOT(8,»0,,-3) 

r 

r RFSFT ORTrJTM 
C 

CALL PL0T<n,,-ll,,-3) 

CALL PLOT(.«i,.«i,-'^) 

r 

r MSP OF L&XIS 

r 

CALL LAXIS(.5,,«»,3,l?HLOG X - AX IS* 12 .6. ♦-2 • 2 ) 

r 

r CRFATF ARRAY TO RF SCALED 
C 

DO 1 1=1, in 

XARRAY{n=FLOAT{ I) 
1 CONTTNUF 
r 

r SCALE DATA USING LOGSCA 
C 

CALL L0GSCArXARRAY,ln,6.,KVAL»NCYCLF) 

r 

r PLOT DATA 

r 

CALL LAXIS(.5»5.,NCYCLF,12HLOG X - AX IS, 12,6. ♦KVAL,2 ) 

r 

r END PLOT 

r 

CALL PLOT(in.,0,,-3) 

STOP 
ENO 
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IDENTIFICATION: AXIS 



S30200 



Purpose : 



This routine draws any length line at 
any angle, divides it into one -inch 
intervals, annotates the divisions with 
appropriate scale values, and labels 
the axis with a centered title. Axis 
is called separately for each axis. 



Restrictions: 



None 



Usage 
Calling Sequence: 



Call AXIS (BGINX, BGINY, IHOL, ±ICHAR, 
AX, AN, FV, DV) 

BIGNX, BIGNY - Coordinates o£ starting 
point of axis line. 

IHOL - Annotation on axis may be an 
alphameric array or hollerth literal . 

+ICHAR - Number of characters in the 
title, an by the sign, determines on 
\Aich side of the line the tick marks 
are made and the title will be placed. 

Note: If positive, all amotation appears 
on the positive side of the axis (y-axis) , 
If negative, all annotation appears on 
the negative side of the axis (x-axis) . 
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Storage Required: 



Method: 



Original Developer: * 



AX - Length of axis line in inches. 
AN - The angle at which the axis is 
drawn. Usually 0" for the x-axis and 
90.0° for the y-axis. 

FV - Starting value to be placed on 
first tick mark. 

DV - Ratio of data units to be added 
per inch of axis. 

This routine is one of the standard 
CALCCM' routines and can be found in 
the Programming CALCOMP PEN PLOTTERS 
manual from CALCCKP, Inc. 

CALCCM', Inc. 
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Identification; Subroutine AXIS2 S30201 

Purpose : 

To draw an axis with control oyerj the number of 
tic marks to be drawn on the axis, the spacing between 
these tic marks, and the number of decimal places on the 
values at the tic marks. 

Usage; F ORTRAN 
Call ing Sequence ; 

CALL AXIS2 CXO, YO, AMAX, A14IN, DELX, AINCH, BCD, 
NCH, NDEC, DELNJ 



where , 



XO, YO Coordinates of the starting point for 
the axis 

AMAX, Maximvun and minimum values on. the' axis 
AMIN (AMIN is the first value, J\MAX the last) 

DELX Increment along axis (the value between 
tic marks) . DELX should divide evenly 
into AMAX-AMIN 

AINCH Length of the axis in inches. If AINCH 
is negative a horizontal axis is drawn 

BCD Hollerith string containing the title 
for the axis 

NCH Number of characters- in the title 

NDEC Nximber of decimals to be used for values 
on tic marks. If NDEC is set to -1, no 
decimals or decimal point is drawn 

DELN 'Value returned by AXIS2 . It contains the 
number of units per. inch for this axis. 
DELN is to be used for scaling data. 
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Connnents 

The user must scale his data to the axis. The 
values of AMAX and AMIN need not be the maximum and 
minimum of the data array. For example, if the data 
ranges from 1000 to 20000^ and the integers 1 thru 20 
are drawn, the values along the axis represent thousands 
and the scaling is: 

XCI) = CX(I) - AMIN*1000)/CDELN*1000) 
Other Subroutines Required: PLOT, SYMBOL, NUMBER 
Storage Required: DECsystem-10 (427) 8 Locations 

Original Source: Tom Talbot/KIIL 

In~house Submitter: Tom Talbot/KHL 

In-house Contact: Tom Talbot/KHL 
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c 

r A TFIT PROfiPAM TO ^HOW THF II55F OF StlBROlfTTNF AXIS2 
r 

CALL PLOTS(n,0,16) 

r 

r .<;ft origin 

c 

CALL PL0T(n,,-Tl.,-'?) 
CALL PL0T(.«i».'5,-3) 

c 

r DRAW X-AXI^»5 TN.LONG^BFGINNING AT COORD. (.5,.S) 

r WITH ? DIGITS TO THF WIGHT OF THF DFCIMAL POINT AT THF TIC MARKS 

r AND WITH THF ANNOTATION ♦X -AXIS' 

r 

CALL AXIS?(,5,,5,5,n,,l,,-6,,8HX - AX IS,fi ,? ^XSC) 

r 

C DRAW Y-AXIS,4 IN,LONG,BFGINNTNG AT COORD. (.5,1.) 

C WITH INTFGEP VALUES AT THF TIC MARKS AND ANNOTATED WITH 'Y-AXIS* 

r 

CALL AXIS2(.5»1. ,2. »-?.♦!., 4. ♦8HY - AXIS»8»-1 tYSC) 

r 

r FND PLOT 

r 

CALL PLOT(R.,n,,-^) 

STOP 

FND 



Identification; Subroutine AXI5 3 S30202 

Purpose : 

To- draw an axis with tic marks at specified 
intervals/ annotate the value of the variable at the tic 
marks and provide an axis identification label. 

Usage; FORTRAN Calling Sequence 

CALL AXIS3 (X,y, THETA, DIST, ORIGIN, DV, TMAJ, 
TMIN, BCD, HGT, N) 



where , 



X,y are the coordinates in floating point 

inches of the starting point of the axis 
with reference to the plotting area origin 
as established by CALPLT. 

THETA is the angle of rotation measured counter- 
clockwise from the X-axis in floating point 
degrees. NOTE: Normally, THETA is 0* for 
an X-axis and 90° for a Y-axis. 

DIST is the length of the axis in floating point 
inches. Should be a multiple of TMAJ. 

ORIGIN is the functional value to be assigned to 
the origin (i.e., the value of the fitst 
scale) , in floating point. 

DV is the adjusted scale factor for the array 
to be plotted (change in value per inch) . 
NOTE: Values of ORIGIN and DV which v/ill 
produce a reasonable scale may be calculated 
using the subroutine ASCALE. 

TMAJ is the distance in floating point inches 
for major tic marks tO.25 inches high). 
Numbers are placed on the axis at the major 
tic marks in accordance with the values of 
ORIGIN and DV. The nvonbers written along 
the axis are adjusted to be between 1000.00 



2.17.2 



TWIN 



BCD 



HGT 



N 



and 0.01 in magnitude. Immediately 
after the last number on the axis is 
placed the caption, xlO^^P, where exp 
is the required exponent. If the 
values are integer multiples, the 
decimal point and decimal places are 
eliminated. A negative TMAJ will cause 
the actual value to be written instead 
of the adjusted value. 

is the divisions per inch in floating 
point for minor tic marks (0.125 inches 
high) . To eliminate minor tic marks 
the following may bfe used: 

TMIN =0. 

is the character Icibel for the axis. 
Can use array . 

is the height of the full-size charac- 
ters in the BCD title. Numbers at the 
tic marks will be (0.75 *I1GT) high. 
HGT is in floating point inches. 

If HGT =0., all annotation will be 
eliminated. 

is an integer specifying the number of 

characters in BCD title. A negative N 

places the annotation on the clockwise 

side of the axis. N = is not allowed. 

If it is desired to have no label, then 
the BCD parameter should be IH, and N 
either + or -1. 



Other Subroutines Required: 
Storage Required: 
Comments : 
Original Source: 
In-house Submitter: 
In'"house Contact: 



PLOT, SYMBOL, NUMBER 
DECsystem-10 (1207) g locations 
Only perpendicular axes 
Langley Research Center 
Tom Talbot 
Tom Talbot 
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r 

r A TEST PROGRAM TO <;HOW THF IISF OF StlBROlJTINE AXIS3 

r 

CALL PL0TS(n,0,i6) 

r 

r <;BT ORTr;TN 

r 

<*ALL PLOT(n,,-n .,-^) 

CALL PLOT{.«>,,5,-3) 

r 

r DRAW X-AXl<;(THrTA = 0,DFr.),BFGINNrNG AT COORD, (O.,0,> 
C S TN.LONG^WITH MAJOR TIC MARK?; FVFPY TNCH AND MINOR TTC 
r MARK:«;,in PPO TNTH. ANNOtatF ON CLOCrWTcF <;rnF OF AXIS 

r WITH ly - AXTSt 

c 

CAl L AXlS?(f^. ♦^•♦0.,5.,n,,l,,l.,10.,flHX - AXlS»0. 14«-8 ) 

r 

C DRAW Y-AXlS(THFTAa9n.DEG) , BEGINNING AT COORD. (0.»95) 

C 6 IN.LONGfWiTH MAJOR TIC MARKS FVFPY HALF INCH AND NO 

C MINOR TIC MARKS. ANNOTATE ON COIJNTERCLOCKWISF SIDE OF 

r AXIS WITH»Y - AXTSt 

c 

CALL AXlST(n, ♦,5,90,,6.,n,,l.,,5,o,,ftHY / AXlS,0»l4»8) 

C 

<" FMD PLOT 

r 

CALL PLOTf 8,,0, ,-i^) 

S TOP 

FND 
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Identification; Siabroutine LAXIS S30203 

Purpose ; 

To draw a log (base 10 axis for use in log-log or semi-log 
plots. 

Usage; FORTRAN Calling Sequence 

CALL LAXIS (XO, YO, NCYCLE, BCD, NCHAR, ALONG, KVAL, IWAY) 

where , 

XO, YO Starting point in inches for axis 

NCYCLE Number of log cycles to be drawn 

BCD Array containing characters for title. 
(H format accepted) 

NCHAR Number of characters in the title. If NCHAR 

is negative, labeling is above axis for x-axis 
or to the right of axis for y-axis 

ALONG Length of the axis in inches 

KVAL The starting exponent of the axis 

IWAY =1 vertical axis is drawn (y-axis) 

=2 horizontal axis is drawn (x-axis) 

Comments : 



The variables NCYCLE and KVAL may be obtained from the routine 
LOGSCA. The three variables NCYCLE, ALONG, AND KVAL must 
agree with the same variables in the routines LOGSCA or 
LSCALE for correct scaling of data. 
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other Subroutines Required: 
Storage Required: 
Original Source: 
In-house Submitter: 
In-house contact: 



PLOT , SYMBOL , NUMBER 
DECsystem-rO (1100) g Locations 
Tom Talbot/KHL 
Tom Talbot/ICUL 
Tom Talbot/KHL 
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A TFST PPOGPAM TO <;H0W THF OSF OF SIJRPOtlT TNF LAXTS 
TALL PLOTSrA,0,"|ft^ 

r 

r <FT ORIfiTN 

r 

CALL PLOTfOf-n,,--?) 
CALL PLOT(«'5,.5,-3) 

r 

r DRAW L0rn-XAXrS»6 T N .LONG, BEr; INNING AT COORD. U5».5) 

C WITH FACH LOG CYCLE 2 IN, LONG» AND STARTING EXPONENT 

C -2. ANNOTATF WITH 'LOG X - AXIS* 

r 



CALL LAXTS(,'>,.S,3,1?HL0G X - AXI S , T? ♦S, »-? , 2 ) 



r 

r DRAW L0n-YAXTS,7 TN,LONr,RFGINNTNG AT COORD, (,5»1,) 

r WITH FACH LOG CYCLF 1 IN, LONGHAND STARTING EXPONENT 

r -->., ANNOTATE WITH 'LOG Y - AXIS* 

r 

CALL LAXIS(,5.1,,7,1?HL0G Y - AXI S, 12 ♦?. »-3 ,1) 

r 

r END PLOT 

r 

CALL PLOT(fl,»n,,-.?) 

STOP 

FMD 
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Identification: Subroutine AXL-ABL S 30204 

<■■■- " ■ ^ ■ I !■ > > ■ I. I II I ■ III ■ ■»■■■■ Il»l. I — ■■■■■■■■I ■■■■■I. '—^1— P^ II ■■ 

Purpose : 

To' number and label the cc-ordinato axes, 
mainly used in conjunction with subroutine GRID. 

Usage: FORTRAN Calling Sequence 

CALL AXLABL (XMAX, YMAX, H, XF , XX, XL, MODX , 

IFMX, YF, YY, YL, MODY, IFMY, IXLABL, 

IX, lYLTvBL, lY) 



where , 



The parameters of the calling sequence are: 

XMAX length of tiie x axis, inches 

Y14AX length of the y axis, inches 

H height of letters a3id numbers, inclj.es 

XF first number to be printed on x~axis 

XX the increment or factor by which successive 
numbers of the x axis differ (positive 
number) 

XL the last nvimber to be printed on the x axis 

MObX the "mode" desired. If I40DX » 1, XX will be 
regarded as an increment (or decrement) , 
resulting in ariDimetic piogression. 
If MODX « 2, XX will be a factor, resulting 
in geometric progression. For example, 
. if XF « 1 and XX = 2, with MODX 
- 1 the successive numijers will be 1, 3, • 
5.... With MODX « 2 they will be 1, 2, 

IFMX a parameter that determines how many decimal 
places will be printed. With IFMX == 0, each 
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number will be printed v/it.k an many decimal 
places as required by the significant 
digits of the number, and no more. With 
IFMX = 1, all the numbers will be printed 
with the same number of decimal places, 
determined by the number that has the most 
places. For example, if XF - 0. and XX 
= .25, with IFMX - the successive numbers 
will be printed 0, 0.25, 0.5, 0.75, 1, 
1.25, etc. V7ith IFMX =• 1 tliey will be 
printed 0.00, 0.25, 0.50, 0.75, 1.00 etc. 

YF analogous to XF, but for y axis 

YY analogous to XX, but for y axis 

YL analogous to XL, but for y axis 

MODY analogous to MODX, but for y axis 

IFMY analogous to IFMX, but for y axis 

IXLABL hollerith matcr.-i.al for labeling the x axi.s. 
I (for example, 16HRANGE, NAUT. MI.) 

IX the nuiriber of hollerith characters in 
IXLABL (in example , IX = 16) . 

lYLABL hollerith material for y axis 

I lY number of characters in lYLABL. 

It was mentioned above that for modes 1 and 2 the 

parameter XX is alv/ays positive. However, XF and XL may 

be either positive or negative (except that in Mode 2 both 

must have the same sign). In Mode 1, if XL > XF, the 

successive numbers will be incremented by the amount XX; 

if XL < XF, they will be decremented by 'XX. 

The maximum number of characters that any one 

number can have is seven (including negative sign, if 
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present, and iucluding deciraaj. point, i£ present) . 
Therefore, if the program call^:; for printing the number 
-5.73928, the number actually printed will be -5.7393. 
The integer 9999999 is the largest number that can be 
printed. The numbers to be prirxted arc alV7ays given in the 
program as floating-point numboro, but V7ill be printed as 
integers (v/itJiout decimal point) if they are v/holc numbers 
and if IFMX (or IFMY) is zero. 

Provision must be made for printing tlie numbers 
and labels below and to the left of the graph by cj-.tabli.i.;hing 
the origin of coordinates above and to the right of the 
lov.'<'r left corner of the plottijig area. For example, 
before doing any plotting, or mailing calls to GRUD, a call 
such as 

CALL PLOT (1.5, 1.5, -3) 
should be made. This provides a 1.5-inch margin in 
which the lettering and nun-ibering will be done. The actual 
size of the necessary margin V7ill of course depend oh the 
size of the letters. 

The user of t:he suyjroutine shou.ld choore a nine (11) 
of numbers .-.uitable lor the total si%e of tlie graph: 
7\ typical appropriate value for II is cxbout 0.02 times XMAX 
or times YilAX, whiciiever is larger. Care mu.st be tiOcen 
that the numbers of the x axis are of si:^e and spacing so 
th;it they do not overleip. 
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Tlic two L'cisic modes permit the axis immberinri 
to go cither in arithmetic or geometric pro<j'ression. A 
third mode (HODX or MODY > 2) provides for ]..lacing 
arbitrary numbers at arbitrary positions on I lie axcr. . 
For example, if MODX =5, the parameter XF becomes on 
array of dimension 5. (It must be so dimenijionod in the 
plotting program - i.e., in this example tlierc must be a 
statement DIMENSION XF(5).) This is the array of numbers 
to be on the x axis. Before calling AXLABL, the program 
naist define the array XF , either by Fortran stcitementi. , 
a RETvD statement, or a DATA stc'.tement .. 

Also^ the parameter XX becomes an array of the 
same dimension. This array defines the positions (in 
inches from the origin) at which thr; numbers of the XF 
array will be printed. 

This mode ca.n be exorcised on either the x or the 
y axis alone, or on both. 

The option that permits one or botli axes to be 
labeled V7ith "powers of ten," instead of with simple numbers 
is obtained by means of a COMMON declaration: 

COMMON/PTEN/IXTEN , lYTEN ' , 
If this declaration is made in the calling routine, and if 
IXTEN is set equal to 1 prior to calling hXLABh, and if 
for example XF =-■ 5, XX =-^1, and MODX •■= 1, then the axis will 
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he labfiDocJ lo'', 10 , .10'' , . . .etc . (instead of 5 , f. , 7 . . . ) . 
Sotting lYTEN -- 1 produces the?, same ro£5ult IJox the y axis;. 
This option is of courEe inost useful for labeling logaritlinic 
scales. 

Tha num]joring and/or descriptive labeling of one 
of the axes can Le omitted, if thin is desired. To oinit 
X-axis numbering, set XX = -- 1. To omi.t y-axir> numbering 
set YY =^ - 1. When this is done, dujvrniy variables must 
be inserted for unused parameters in the Ccilling statem:-nt 
e.g., for XF , XL , and MODX vjhen XX = - 1 . 

If it is desired, for some reason, to omit the 
x--iixis descriptive .label, tl'iis result is ciccoraplished by 
setting DIX = 0. Hov/ever, a dummy parameter must of course 
then be inserted for IXLABL — e.g., IIIX, If it is dcsJrcd 
to omit the y-axis descriptive label, this is done by setting 
lY -- 0, £ollov;ing a dummy parameter lYLABL. 

Subroutine AXLABL has subordinate subxoutines 
named CErvia'KR. DECDGT, ixnd PV'RTI:;N, 



Other Subroutines Required; 

Storage Required: 
Original Source 



PLOT, SYMBOL, IHUIBER, CLliNiTJJR, 
DECDGT, PVffiTI^N 



Computer Subroutines to plot 
and label curves for Scientific 
Naval Research Lab., Washington 
0,C. 
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In-house Contact: 



Tom 'i'ciji/oL/iaiJi 
Tom Talbot/KIIL 
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r t TF«^T PPOnPAM TO qHOW THF \)<^P OF «;i)RPnUTTNF5; GRTD AND AXLARL 

DIMFN«;iON XX(5) ,XF(5),YYf 1 ) fYPd) 

COMMON « T PFN M XTFN ♦ I YTFN 

r 

CALL DLOT^fn,n»Tft) 

r 

r cFT oQTriTN 

r 

CALL PL0T(n.,-ll,»-3) 
CALL PLOT(4.,l.,-->) 

r 

r LAYOUT «»x'=; rPTn,MO TTr mapk<; 

r 

CALL GRID(5.,5.,n. ,n.,5,5,l,l,3,3»l,l,l,l) 

r <;ft papamftfp*^ 

r 

xFn)=n, 
xxfn = i. 

YF(1 V=?* 
YYM 1 = 1, 

r 

r M^F «;ttRPOIfTTNP AXtAPI. 

r 

CALL AXLABL(5.,'S,,.14,XF,XX,5.,1,1,YF,YY,7.,1,0,8HX - AXT.S,8,8HY 

* AXTCR) 

r 

r FND PLOT 
r 

CALL PLOT(R, jO,,--:*) 

r 

r CFT NFw r;PTn 

r 

CALL CRTOfA. ♦5,,. OR, ,n8,10, 2,1, ?♦?,?, 1,1, 1,10) 

r 

r FTLL ARPAV FOR ANNOTATTOW 

r 

no 1 T=r,s 

XF(I)=FL0AT(I)-3. 
XX(T) = FLOATn)-l. 
1 CONTTMHF 
TYTFMsl 

r 

CALL AXl ARL(4.,«^,,.14,XF,XX,4.,5,n,YF,YY,4.,l,l,fiHX - AXl.S,8,BHy 

* AXT<;,P) 

r 

r FNH PLOT 
r 

CALL PL0T(R,,n,,-3) 
C 
r CFT MFW r,PTn 

CALL GRID (4,, 5,,-. n ft, -.08 ,10, 2,1,?,?, 2, 1,1 ♦1,10) 
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r ■ ■ 

r rrr D/vDAMrTFPC 

yyn) = ?, 

vYn) = i, 

TYTFM=] 

r ■ 

CALL AXLABL(4.,5.,,14,XF,XX»5.,2,0,YF,YY,4,,1,1,8HX - AXl<;»8,flHY - 
♦ AXT^,R) 

r 

r TMH PLOT 

r 

TALL DLi^Tf R.,n,,_^^ 

c;TOP 

FMH 
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Xdentification; Subroutine LOGAX S30205 
Purpose ; 

Plots axes with log^Q scale. Options permit 
labeling on either side of axis and grid lines at powers 

of 10. 

Usage; FORTRAN Calling Sequence 

CALL LOGAX (XO, YO, NCYCLE, BCD, NCHAR, AXLEN , 
KVAL, GRID) 



where , 



XO, YO Coordinates at origin 

NCYCLE No, log cycles along axis 

BCD Array containing characters for title 

NCHAR No. characters in title; positive for 
conventional labeling, negative for 
opposite side of axis. 

AXLEN Axis length, inches; positive for Y-cixis 

(vertical), negative for X-axis (horizontal) 

KV/iL Pov/er of 10 at origin 

GRID Length of grid lines (use length of other 
axis); if zero, conventional tics are 
used - no grid. 



Other Subroutines Required: 
Storage Required: 
Conmients : 



Restrictions : 

Original Source:' 
In~house Submitter: 

In-house Contact : 



PLOT, SYMBOL, NUMBER 
DECsystem~10 cere requirements 

Modification of LAXIS (Improved 
to adjust axis labeling according 
to length of one cycle and provide 
grid option) 

Do not use less than 1 inch per 
cycle 

Alex Robb, Tom Talbot/KHL 
Alex Robb/KUL 

Alex Robb/KHI. 



CALL PL0T8<0r0»lA) 

CALL PLOT (0. ,-11., -3) 

CALL P^OT (5. . .t;,-3> 
i WRIT£0,29) 
29 FORMAT (' ?') 

READ(5»130,E:Nn = 5l?.) LOC , NCYCLE»KVAL , X0» Y0, XLEN, XTIC 

IF(LOC.GT.I) GO TO 53 
CALL LOGAX(X0,Y0,NCYCLE,8HX-AXIS , 8»L0C# XLENi KVAL, XT IC) 

WRITE(5,29) 
■READ (5. 120, END =50) LOC . NCYGLE ,KVAL , Xfl. Y0. YLEN, VTlC 
CALL lCGAX(X0,Y0.NCYClE,8HY - AXI S , 8«L0C. YLEN» KVAL. YT IC) 
CALL PlOT(ABS{XLEN)*3. .0.1-3) 
GO TO 1 
50 CONTINUE 

CALL PLOT (10. ,0. ,-3) 
STOP 
lt»0 F0RMaT(3I»4F) 

END ■ , . 
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where , 



T.dent.i.J:ication; gubroutine DIMEN S S 3 2 6 

Purpose ; 

To. draw dimensions for drafting purposes 

Usage; FORTRAN Calling Sequence 

CALL DIMENS (XO, YO , DS, THETAt SCALE) 



XO, YO are the floating point coordinates of the 
left end of the dimension. 



DS 



is the floating point length of the 
dimension. 



THETA is the angle in floating point degrees at 
which the diiiTension is to be drawn. 

SCALE is the scale of the drawing in floating point. 



Other Subroutines Required: 
Storage Required : 
Restrictions: 



PLOT; NUMBER 

DECsystem~10 (406) „ Locations 

The dimensions are limited to 
three (3) types which are 
determined automatically by 
the actual plotted length of 
the dimension. Type 1 is 
utilized when (DS*SCALE) 
is <_ 0.80 inch, Type 2 
when 0.81 <_ 1.19 inches, and 
Type 3 when 1.20 inches or 
greater . 



Type l-t> 
Type 2|^|- 
Type 3|^ 
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Original Source: 
In-house Submitter; 
In-house Contact: 



Langley Research Center 
Tom Talbot/KHL 
Tom Talbot/KHL 
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A TFST PROGRAM TO ,^HOW THF USF OF SUBROIITTNES DIMENS AND LIGRID 
CALL PLOT«;(n,o,i6) 



r 

r CFT ODTGTM 

r 



<-ALL PLOT(r^.,_^l,,_-,) 
Ch\ \ Pf.nT(,^,,B;^_-, , 



DRAW LTNFAR GRTD LTNF^ VFPTTCAL L INF5 APF ^ INCH APART(XS=1.) 

rAM . Tr-r,T,^,« « HORIZONTAL LINPS ARF 2 INCHFS APART(Y^ = 2.) 
CALL Lrr,RID(n,,0,,l.,2.,5,3) 

DRAW 3 TYPF«^ OF RFPRF<^FNTAT I VF niMFN^ION.c 



CALL niMFN«^n..6.?,l,,n.,i,) 
CALL nTMCMC(S.?,?.,?.,on.,i,) 
CALL PLOT(4,'i,ft,,3^ 
CALL PLOTf4.'?,'^,7«^,?) 
CALL DTMFN<;(A.,6.2,,5,n,,i,) 
r 

r PND PLOT 

f" CALL PLOT( ft., 0, ,-■:() 

<^TOP 

FMn 
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Identification; Subroutine PAXIS S3 020 7 

Purpose t 

To draw a probability axis with tic marks at specified 
intervals and annotate the value at the tic marks, when space 

permits. 

Usage ; FORTRAN 

Calling Sequence; 

CALL PAXIS (XO, YO, HIGH) 



where , 



XO, YO Coordinates of the starting point for the 
axis (in inches) 

HIGH Height of the axis in inches 
Comments ; 



The user must scale his data to the probability axis 



drawn , 



Other Subroutines Required; PLOT, NUMBER 

Storage Required; 27 4 g 

Original Source; Scott Moffatt/KHL 

In-house Submitter; Tom Talbot/KHL 

In-house Contact; Tom Talbot/KHL 
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r 

r A TF<;t PP0<;PAM to <;H0W THF tl.^F OF SUBROUTINE PAXIS 

r 

CALL PLOTS(<>in,iA) 

r 

r >FT ORITiIN 
r 

CALL PL0T(n,,-il.,-3) 

CALL PLOT(0.»l,»-3) 

r 

r «FT HFIGHT OF PPORAPILITV AXI.«^ 

r . 

HHsft.rt 

.r 

r <;FT AXIC, rOOPOTNATPC 

r 

XOsA.n 
YO=n,r» 
CALL PLOT(3.,0, ♦-■:?) 

r 

r PLOT PPORABTLTTY AXTc 
C 

CALL PAXT.S(XO,YO,HH) 
r 

r <:FT OPir.TN POP NFXT PLOT 
r 

CALL PLOT(3,,0,»-^y 

r 

HH»4«0 

CALL PAXT«;(XO,Yn,HH) 

CALL PLOTC^.tO,,-^) 
C 

HH=?,0 

CALL PAXI<;(XO,vn,HH) 
C 
r FND PLOT 



r 



CALL PL0T{3.,r>,,-.3) 

«;top 

FND 
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IDENTIFICATION: SUBROUTINE GRIDD S30210 



Purpose : 

To draw grid lines inside specific boundaries with an 
optimal box left open for labels. GRIDD can be used with 
AXIS routines to add the finish lines to the outside 
parimeter of a graph. 

Usage; FORTRAN 

Gaining Sequence: 

CALL GRIDD (XO, YO, XGRID, YGRID, KX, KY, LAST, lOPN, 
XLOW, XUP, YLOW, YUP) 



where , 

XO, YO 



The coordinates of lower left hand corner 
of the graph. 



XGRID, YGRID The lengh of the X and Y axis. 

KX, KY The number of cycles on each axis. 

One less line than the number will be 
drawn. 



LAST 



lOPN 



XLOW, YLOW 
XUP, YUP 



Options to draw graph outlines. 

= - No outlines drawn 

= 1 - Right outline drawn 

= 2 - Upper outline drawn 

=3 - Both outlines drawn 

This place the Box option 

= - No box left 

= 1 - Box with no borders 

= 2 - Box with border drawn 

Bottom coordinates of the box boundaries 
relative to XO, YO. 

Top coordinates of the box boundaries 
relative to XO, YO. 



other Subroutine Required; 
Storage Required: 
Restrictions: 

Original Source: 
In-house Submitter: 
In-house Contact: 



Plot 

DECsystem-10 (3060) « Locations 

If IPON = 0, XLOW, YLOW, XUP, 
YUP must contain dummy 
variables . 

Langely Research Center 

Tom Talbot, KHL 

Tom Talbot, KHL 
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Identification: Subroutine GRID S 30 211 

Purpose : 

To draw different combinations of linear and 
logarithmic grid lines with decade, unit and fractional 
tic marks. 

Usage : FORTRAN 

Calling Sequence ; . 

CALL GRID (XMAX, YMAX , TICX, TICY , NCX , NCY , IJ-.X , 
LLY, MODX, MODY, MT , MF, NI , NF) 



where / 



XMJiX The desired width of the coordinate grid, 
inches (X-axis) . 

YMAX The desired height of the coordinate grid, 
inches (Y-axis) . 

TICX The length of the X-axis tic marks, if 
any. (If none are called for, enter 
zero). If TICX is positive, tic marks 
will be directed inward from border; if 
negative, they will be outward* 

TICY The length of the Y-axis tic marks. 

NCX The number of linear decades or log cycles 
in the X direction . 

NCY The number of linear decades or log cycles 
in the Y direction. 

LLX An integer that- indicates whether the X- 
coordinate scale is to be linear or 
logarithmic? 1 for linear, 2 for logarithmic, 
3 for reversed-logaritlimic. 

» 

LLY An integer that indicates whether the 
Y-coordinate scale is to be linear or 
logarithmic; 1 for linear, 2 for logarithmic, 
, 3 f or reversed-logarithmic. 

MODX An integer indicating whether or not 
tic marJrs are desired and v/hether or 

not fine-pea coordiaate liueii arc desired 
for the units in each decade, in the X 



5-17.11 



dimension of the grid. (A lifsting and 
de script ion of the various modes is givoi 
on the next page . ) 

MODY Parameter corresponding to ,i<50DX for the 
Y dimension of tne grid. 

MI An integer between 1 and 9 that indicates 
at which unit the initial decade of an 
X-axis logaritluaic scale ia to begin. 

MF An integer between 2 and 10 that indicates 
at vj-hich unit the final decade of an X-axis 
logarithmic scale is to end . 

NI Same as Ml except refers to' Y-axis. 

NF Same as MF except refers to Y-axis. 
(The last four parameters are not used by the subroutine 
when a linear-scale grid is plotted, but "dummy" parameters 
must be inserted in the parameter Id.st in this case.) 

With respect, to logarithmic decades, the terms 
decade, unit, and fraction are defined as indicated by 
the following diagram: 



dec«idG 

r—. '- ^ 

• unit tick marks * • 

Nfraotjonal ticks 

Lu..LU.iiii.ljiji :iiiili!fitiiii!H!i;!iii[|j|iliiiilii!iiiiM I I 1 1 nniliim niiliiiiniiil ipiiil i|iilii|ii.. 
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Modes 

The modos (ivlODX, MODY) that can be called arc 

listed below. The nximerals in parenthesis refer to 

Gerber penturret locations: 

i4ode = 1: (For log scales only) Decade lines, 
broad pen (12) . 
Unit lines, fine pen (10). 
Fractional tic marks, fine pen (10). 



Mode = 11: 



Mode = 2 



Mode - 3 



Mode = 4 



Mode .= .41', 



Mode " 5; 



(For log scales only) Decade lines, 
broad pen (12) 

Unit tic marks, fine pen (10). 
No fractional tic marks. 

Decade lines, broad pen (12). 
No unit lines. 

Fractional tic marks, fine pen (10) . 
For log scales, tic marks at unit 
intervals are 2.0 times as long. 
For linear scales, there are 10 tic 
marks per decade, vrith fifth tic 
mark 2.0 times as long. 

Decade lines, broad pen (12) . 

Log scales, unit lines fine pen (10). 

Linear scales, 10 fine lines per 

decade (10) . 

Log scales, no fractional lines. 

No tic marks. 

(For linear scales only) Decade lines, 
broad pen (12) 

Four tic marks per decade (5 intervals), 
fine pen (10) . 

(For linear scales only) Decade lines, 

broad pen (12) . 

Four fine lines per docade (10) . 

No tic marks. 

(For log scales only) Decade lines, 
broad pen (12) . 
Unit lines, medivim pen (11). 
Fractional lines, fine pen (10). 
No tic marks. 



Comineiits 

Grid Opening Option 

An opening can be left in the grid lines to 
annotate the curves or for any other purpose. 

The following common statement should be included 
in the calling routine. 

COMMON/OJ?N/XA , XB, YA, YB, lOPN 

Prior to calling Subroutine GRID, the statement 

lOPN = 1 is made, XA, XB, YA and YB are given numerical 

values (inches) . 

XA X coordinate of left edge of the grid 
opening 

XB X coordinate of right edge of tlie grid 
opening 

YA . Y coordinate of the bottom edge of the 
opening 

YB Y coordinate of the top edge of the opening 

The opening can be surrounded by solid lines 
(rectangle) if desired or it can be left "OPEN". 

If the lines are desired, a call to subroutine 
BOX is made (CALL BOX) just after the call to GRID. 

If, after a call to GRID with the "opening" 
option, another call v;ithout an opening is' desired, sot 
lOPN = to deactivate the opening feature. 

Other Subroutines Required: PLOT, BOX, COORDS, TICKS, 

PENCHG 

Storage Required: " DECsystem-10 

(301^0) (. locations 

I 



'I n .^ 



Restri-ctions 



Original Source: 



In --house Submitter 
In-house Contact: 



The user must provide a dummy 
routine 

e.g. Subroutine PENCHG 

SUBROUTINE PENCHG (IPEK) 

RETURN 

END 

Computer Subroutine to Plot and 
Label Curves for Scientific 
Publications , 
Naval Research Lab. 
Washington, D.C. 

Tom Talbot/KHL 

Tom Talbot/KHL 
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r A tf5;t ppogpam to <;Hnw thf n«;F of «;tipPOUTTNF grid 

r 

rOMMOM/OPN/y? ,y?,Yi ,v?,topm 
r • ,^FT COORD. FOR GRTD OPFN ? NG ♦ I OPN= 1 ACTrVATF.S THF GRID OPFNING 
DATA Xl,X2»Yl,Y2,rOPM/3.,4.,3.,4.,l/ 
CALL PL0T<^(0, 0,1ft) 

r 

r <.V1 ORIGIN 

r 

CALL PL0T(0,,-n,,_3) 

CALL PLOTf^,,!,,--^) 
C 

r IIc;F GRin ,<;t(RROltTTNr 

c 



c 

c 



CALL GRID(5,,5.,0,,n.,5,5,],l,3,3,l,]«l,l) 



r 

C CALL TO .SUBROUTINF BOX DRAWS LINES AROUND T.HE OPENING 

CALL pox 

r 

C CALL TO STANDARD ROUTINE SYMBOL WILL ANNOTATE THE BOX 

r 

CALL «^YMROL(3,?S,3.?'5,.2S,?H ♦n,,3) 

r 

r FND PtOT 

c 

CALL PLOT(8,»n,,-3) 
C 

r MAtfF NFW GRID 

<" ! 

CALL GRID^5.,5.,.n8,,n8,l,l,2,2,l ♦1»1»10,1»10) 
CALL PL0Ttfl,»0,,-3) 
r IOPN»0 DF-ACTTVATF*^ THF GRID OPFNT^'G 

CALL GRID(5,»5,,-,08,-.nR,i,i,2,2»l,l,l,10,l,lb) 
CALL PLOT(8.,n,,-3) 

lopN=n 

CALL GRiD(4.,'S,,-,ne,_,nft,in,2»l,?»?,2,i ,1,1 ,10) 
CALL PLOT(8.,0,,-3) 

c 

CALL CjRlD(5,»4«,0.,n,,l,l,2,2,3,3»l,in»l,10) 
CALL PLOT(8,.n,,-3) 

r 

CALL GRID(6»»5,,,08,.n8,5,'^,l ,l»A»41,l,],l,n 
CALL PLOT(ft,,0, ,-■:)) 

C 

CALL GRTn(6,,«>,,-,02,-,o?,9,l ,2,2, «5,«5,"; ♦3,2,0 j 

C 

r PND PLOT 

C 

CALL PLOT(8.»n.,-3) 

STOP 

END 
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area, 



Identification; Subroutine LIGRID S30212 

Purpose ; 

To draw linear grid over the specified plotting 



where , 



Usage; FORTRAN Calling Sequence 
CALL LIGRID (X, Y, XS, YS , M, N) 



X, Y 



XS 



YS 

M 

W 
Other Subroutines -Required 
Storage Required: 
Original Source: 
In-house Submitter: 
In -house Contact: 



are the page coordinates in floating point 
inches of the origin (lower left-hand corner) 
of the grid. 

is the increment in floating point iiiches 
between vertical grid lines. 

is the increment in floating i:)oint inches 
between horizontal grid lines. I 

is the number of increments along the X-axis. 

is the number of increments along the Y-axis. 

PLOT 

DECsystem-10 (153) g Locations 

Langley Research Center 

Tom Talbot/KHL 

Tom Talbot/Kim 



2.17.12 . 



SPECIAL SYMBOLS 
S303XX 



2.18 



IDENTIFICATION: 



SYMBOL * S30300 



Purpose: 



Restrictions : 



This routine produces plot annotation 
at any angle and in practically any 
size. There are tvro SYMBOL call for- 
mats: (1) the standard call vfliich can 
be used to draw text; and (2) the 
special call v^ich can be used to draw 
special symbols. 

None 



Usage : 
Calling Sequence: 



Call SYMBOL (BGINX, BGINY, HI, IHOL, 
AN, TNCH) 

BGINX, - Coordinates, in inches, of the 
lower . 

BGINY - Left comer of the first character. 

HT - Height, in inches, of character to 
be plotted (e.g., .07, .14, .21, ...) 

IHOL - Text to be used as annotation. 
Usually stored in BCD or A-type format. 

AN - "Hie angle in degrees, at vrtiich the 
annotation is drawn. 



2.18.0 



INCH - Number of characters to be plotted 
from IHOL. IF NCH > 0, data must be 
left-justified in IHOL. 

IF NCH = 0, one alphanumeric character' 
is produced. 



Usage: 
Calling Sequence: 



Storage Required: 
Method: 



Call SYMBOL (BGINX, BGINY, HT, INT, AN, 
-IC) 

BEGINX - Coordinates, in inches, of the 
lower left. 

BGINY - comer of the first character. 

HT - Height, in inches, of the character 
to be plotted. 

INT - Integer equival«it of the desired 

symbol. 

AN - The angle, in degrees, at vdiich the 
annotation is dravn. 

-IC - Determines the per status. 

- -1, Pen up during move, symbol produced. 

« -2, Pen dovn during move, symbol produced. 

This routine is one of the standard CAIXXM' 
routines and can be found in the mOGRAMONG 
CALCCMP PEN PLOTTERS manual from CALCCM>, Inc. 
2.18.0.1 



} I 



Original Developer: * CALCCM>, Inc. 



2.18.0.2 



IDENTIFICATICN: 



NUMBER 



S30301 



Purpose: 



Restrictions : 

Usage: 
Calling Sequence: 



Note: 



This routine functions as a pre-processor 
to the symbol subroutine. It converts a 
real variable to the appropriate fixed- 
point equivalent so that it may be 
plotted. 

ND should be less than or equal to nine. 

Call Number (BGINX, BGINY, HT, FN, AN, 
IND) 

BGINX, Coordinates in inches, of the 
lower left comer BGINY of the first 
character. , 

HT - Height in inches of character to 
be plotted (.e.g., .07, .14, .21,...) 

FN - Location of the floting-point num- 
ber that is to be converted and plotted. 

AN - The angle in degrees, at vdiich 
the number is drawn. 

±ND - Determines the prevision of the 
number FN. 

If 

ND > = > the number of digits to the 
right of the decimal point. 



I / 



ND = = > plots, rounded integer and 
decimal point. 

ND < = > digits are truncated after 
rounding. 



Storage Required: 
Method: 



This routine is one of the standard 
CALCC3MP routines and can be found in 
the PROGRAMING CALCOMP PEN PLOTTERS 
Manual from CALCCWP, Inc. 



Original Developer: 



CALCOMP, Inc. 



2.18.1.1 



Identification ; Subroutine ENUMBR 30302 

Purpose ; 

To convert a floating point nxanber to ASCII (expressed in 
scientific notation) , and draw the resulting alphanumeric 
characters. 

Usage ; FORTRAN Calling Sequence 

CALL ENUMBR (X, Y, SIZE, FPN, THETA, N) 



where , 

X,Y 



are the coordinates in floating point inches 
of the left lower corner of the first digit 
of output. 

SIZE is the height of the plotted number is floating 
point inches (see symbol routine) , 

FPN is the floating point number to be drawn. 

THETA is the angle in floating point degrees at which 
the number is to be drawn (see symbol routine) . 

N is an integer specifying the number of significant 
digits to be drawn. " 

If N is negative, the digits will be to the 
right of the decimal. 

Examples: 

FPN = 9093.1200 

Ci^LL ENUMBR (X,Y,S ,PPN,0. ,3) 
wxll draw; 

909 X 10^ 



CALL ENUMBR (X,Y,S,PPN,0. ,-3) 
wTII draw; 4 
.909 X 10* 



2.18.2 



FPN - .0.1020 

CALL ENUMBR (X , Y , S , FPN , . , 2 ) 

will draw; 

1 X 10~^ 

CALL EKllMBR iX,Y,S, FPN , . , -2 ) 
will drav:; 



.1 X 10 



_T 



Other Subroutines Required 
Storage Required : 
Restrictions : 

Original Source: 
In -house Submitter: 
In-hoxise Contact: 



SYMBOL, NUMBER, WHERE 

DECsystem-10 (457) g Locations 

The resulting nuinber is 
restricted to a maximum 
of 8 significant digits. 

Langley Research Center 

Tom Talbot/KHL 

Tom Talbot/iaiL 



O TOO T 



r 



C A TF<;T PROfiPAM TO SHOW THF t/SF OF SUBROUTINES F-NUMBR AND FNUM 
CALL PLOTS(n, 0,16) 



r 

r SFT ORTr,TN 

r 



CALL PL0T(o,,-ii,,_3) 
CALL DLOTf.S,. «,♦-::<) 

r 

r PROVIDF THC FLOAT TNG POTNT NUMRFR FPN(FM) 

r, 

FPM = 1?^6«i. 

r 

r TALL SUROOMTINFS 

r 

CALL F NUMBP ( 2 . , 2 . ♦ • 28 , FPN ♦ • , 3 ) 

CALL ENUMBR(2.,4.,.14,FPN,0.,4) 

r 

CALL PNUM(2.,6,,,l<Sf,FN,9n,,5) 

r 

r FND PLOT 
r 

CAI L PLOTf fi,,n,,_^) 

STOP 

FND 



#»L*.l?^ 



t f 



o 
O 



in 

CO 
CS) 



1234 X 10» 



123 X 10^ 



Identification ; Subroutine LABEL S30303 

Purpose ; 

To draw a label following a given curve defined by 



points (X^r Y^) . 



Where, 



Usage ; FORTRAN Calling Sequence 

CALL LABEL (X, Y, NTOT, XM, DELX, YM, DELY, K, 
LAB, NLBL, HT, IWAY) . 



X,Y Arrays defining curve 

NTOT Number of points in arrays X and Y 

XM Minimum value for X-array 

DELX Units per inch in X-direction 

YM Minimum value for Y-array 

DELY Units per inch in Y-direction 

K Index in X-array to start label 

LAB Axray containing label (H format accepted) 

NLBL Number of characters in label. If NLBL 

is negative, the label is drawn below the 

curve . 

HT The height of the characters to be drawn 
IWAY = 0, characters are drawn vertically 

» 1, slope of characters follow slope of curve 
(in both cases, the title will follow the curve) 



2.18.3 



Comments ; 

LABEL draws a label following a given curve 
defined by the arrays X and Y. The label will start at 
X(K), Y(K) and v/ill be slightly above or below the curve. 
The label will not continue beyond the end of the curve 
and the curve is not drawn. 



Subroutine required: PLOT 
Other Subroutines Required: 
Storage Required: 
Restrictions : 

Original Source: 
In-house Submitter : 
In-house Contact: 



PLOT 

DECsystem-10 (311) g Locations 

The label will not continue 
beyond the end of the curve . 

Tom Talbot/KIIL 

Tom Talbot/KIiL i 

I 
Tom Talbot/KHI. i 



2.1^.3.1 



r 

r A TF«;t PROfiPAM to «;H0W THP USF of SIIRPOUTINF LARFL (FLIPS) 

nTMFN<;TON X(100)»V(100) 

r 

r CFT nPT^TM 

r 

CALL PLOT{0.»-n..-^) 
TALL OLOT(,R,,«;,-35 

r ■ 

C DRAW X AND Y AXES 

r 

CALL AXl.<;2(0..0.,3,,n,a.»-4.,8HX - AXI <;,8,-l ^XSC) 
CALL AXl5;2(0.»n,,5.,O,,l.» 6., SHY - AXI,<;,8,-1 ,YSO 

r GFNFRATF ARRAY*; FOP ruRvF<; 

r 

DO 1 T=1 *inn 
X(T)=FLOAT(I)»,0':i 

Y(T)«:7.*<TN(Xr TM 
1 rONTTN'UP 

r 

r DRAW A SINF CURVF 
C 

CALL ALINE{X»Y,inO,0,,xSC,0.,YSC) 

r 

r LARFL THr c;TNF riJRVF 

r 

CALL LARFL(X.Y,inn,n.,xSC,n.,YSC,ll,ir>HY = STN(X ) ,10,, 12 , 1 > 
r cnPFPTMDO.CF AN FLLTPf-F 

r 

CALL FLTPS(?.5,3,,l,,2.,3n.,0,,36f).,'^) 

r 

r FND PLOT 

r 

CALL PLOT(8.,0.,-3) 

•^TOP 

FMO 



91ft •» 9 



I i 




Identification: Subroutine PARROW S30305 
Purpos e; 

To. draw an arrowhead at the end of the line 
proportional to the length of the line. 

Usage; FORTRrvN Ca l ling Sequence 

CALL PARROW (XA, YA, XB, YB, NC, ALUGTIl) 

where, 

XA,YA are the starting coordinates of the 
line in floating point inches. 

XB,YB ■ are the coordinates of the arrow- 
head at the end of the line. 

NC =1 plot the line from A to B and then 

draw an arrowhead one-fourth of 
the length of line AB. 

=2 draw ah arrowhead at B one-fourth 
of the length of line AB; 

=-1 plot the line from A to B and then 
draw an arrowhead as specified 
I by ALNGTII . 

=-2 draw an arrowhead at B of length 
as specified by ALNGTH. 

ALNGTII is the length of the arrowhead whcm 

NC is negative; it may be omitted when 
NC is positive. 

Commen ts : . 

The PARROVJ routine allows the plotting of an 

.arrowhead that is either proportioned to the length of 

line it is associated with or is of a stated length (/vLNGTH) 

The arrowhead usod subtends an angle of. 60®, It does not 

use. the LINE routine. i 



•? 1 9 a 



I 



I 
I 

'■ other Subroutines Required: PLOT 



Storage Required: DECsyste.m-10 t237)g Locations 

Original Source: Langley Research Center 

In-house Submitter: Tom Talbot/KHL 

In-house Contact: Tom Talbot/KllL 



O 1 U /I T 



r ■ ^ , 

r A tf5;t pRor,PAM to <;how thf usf of SDPRonTiNE parrow(L<;<:alf, 

r LAXIC,AXI<;2#<^YMR0L) 

r 

DIMFN«^TON Xnon),Y(inO) 
r&l.L PLOT.«;fn,0,if,) 

r 

r <;ft ORir,TN 

r 

rh\A_ oLOT(n,,_n ,,-.^) 

CALL OLOTfiC;, ,«;,-■:<) 
r 

r CFT OATA AOOAY*; 

r 

no 1 T*l,li^n 
X( I)«FL0AT(I)*.n«5 
Y(n=?rN{X(T) ) 
Y(T)=FXP(Y(T) S 
\ CONTINUF 

r 

c <;ft valuf«; for log axt<; 

r 

(fVAl.s-? 

MrYri.Ps'a 
r 
r <irALP" DATA 

r 

CALL LSCALF(Y,inn,KVAL,4,,NCYCLE) 
r 
r DRAW LOG AXI.^ FOR Y ARRAY 

r 

CALL LAXI5;(0,,r>,,NCYrLF,15HY = FXP ( 5;TN r X ) ) ♦ 1 5 »4. ,KVAL,n 

r , 

r DRAW LTNFAR AXIS FOR X ARRAY 
r 

CALL AXT<:?('^.,'^.,«5,,n,,l.,-6,,lHX,l«-1 ,V<;C) 

r 

r DRAW CURVF 

r 

CALL ALTNF(X»Y,lf>f^»n,,X5^C,'^.«l. ) 
r 

r DRAW ARROW TO ANNOTATF rjjRVF 

C 

CALL PARR0W(3,,4,,2.,3,3,-1,.2) 

CALL .'=iYMR0L(3.1,4,U,14,1«iHY = FXP( r,TN ( X ) WO, , 15 ) 

r 

r FMD Of.OT 

CAI L PLOT(R,,n,,--^) 

<;top 

FMD 



2.18.4.2 



Y = EXP(SIN(X)) 



X 

z 

I— t 

CO 

a. 

X 

bJ 

n 

>- 




10"i 



5+ 
4 



10 



H 

5 



IDENTinCATION; SUBROUTINE ENUM S30306 



Purpose ; 



To convert a floating point ntunber to ASCII (expressed 
in scientific notation) , and draw the resulting alpha- 
numeric characters. 

Usage; F ORTRAN Calling Sequence 

CALL ENUM (X, Y, H, FN, A, N) 

where, 

X,Y are the coordinates in floating point 
inches of the left lower corner of the 
first digit of output. 

H is the height of the plotted number in 
floating point inches. 

FN is the floating point number to be drawn. 

A is the angle in floating point degrees at 
Ti^hich the number is to be drawn. 

N. is an integer specifying the number of 
significant digits to be drawn. 

if N is negative, the digits will be to the 
right of the decimal. 



ExEunples : 



FN = 9093.1200 

CALL ENUM (X, Y, S, FN, 0.,3) 

will draw; 

909 X 10^ 

CALL ENUM (X, Y, S, FN, 0.,-3) 
will draw; 

.909 X 10^ 



FN - .01020 

CALL ENUM (X, Y, S, FN, 0.,2) 

will draw; 

1 X lO"^ 

CALL ENUM (X, Y, S, FN, 0., -2) 
will draw; 

.1 X lO"""" 



Other Subroutines Required; 
Storage Required: 
Restrictions: 

Original Soruce: 
In-house Submitter: 
In-house Contact: 



SYMBOL, NUMBER 

DECsystem-10 (457) „ Locations 

The resulting nvunber is 
restricted to a maximum of 8 
significant digits. 

Langley Research Center 

Tom Talbot 

Tom Talbot, KHL 



2.18.5-1 



\ ! 



SPECIAL FUNCTIONS 



.q304XX 



2.19 



SPECIAL FORMS 
S305XX 



2.20 



p 



Identification; Subroutine ARC S30501 

Purpose; 

To draw a circle or portion there of 

Usage ; FORTRAN Calling Sequence 

CALL ARC {.\^,yO, TgO/ TliF, RO, IP) 




where , 



M,yv 



THO 



THf 



RO 



IP 



are the page coordinates of the center 
of the circle, 

is the initial angle between the radius 
and the X-axis in floating point degrees. 

is the final angle between the radius and 
the X-axis in floating point degrees. 

is the length of radius in floating point 
inches . 



= 2 move to start of arc segment with pen 
down. (Trace a line to starting point). 

- 3 move to start of arc segment with pen up. 

Other Subroutines Required: PLOT 

Storage Required: DECsystem-10 (174) g Locations 



Restrictions : 
Original Source: 
In-house Submitter; 
In-house Contact: 



No plot if THE ■ THO 
Langley Research Center 
Tom Talbot/KHL 
Tom Talbot/KHL 

2.20.1 



r 

r A TFe;T PROGRAM TO ,<^HOW THE USF OF SURROUTINE ARC 

CALL PL0TS(0,n,16) 

r 

r ePT opTr;TN 
r 

CALL PLOT(n,,-n ..--s) 

CALL PL0T(l,»l.,-3j 

r 

r ORAW X ANH Y AXFc; 

c 

CALL AXIfi2(n,,0,,6.,0., + 6.,8HY - AX IS ,8,-1 , YSC) 
CALL AXIS2(n.,n,,3,,o,,-A.,8HX - AX T<; ,8,-l ,X5ia 

C DRAW APr(riRCLF) 

r 

CALL ARC(2.,?.,n,,?6n,,l.5,-?) 

r 

r FNin Pf OT 
C 

CALL .PLOT(B,,0, »-3 ) 

STOP 

FNO 



•» •>/* ^ ^ 



5-r 



to 

X. 

cr 



2-- 



1 -- 








X - AXIS 



2.20.1.2 



\ ■ 



where , 



Identification; Subroutine. BLIPS S30502 

Purpose ; 

To draw an ellipse or portion thereof. 




Usage; FORTRAN Calling Sequence 

CALL ELIPS (XO, YO, A, B, fiLP^A, THETAO, THETAF , IV) 



XO, YO are the page coordinates of the starting 
point on the ellipse in floating point 
inches. This is not the coordinate, of 
the center or of a fucus of the ellipse. 

A is the length of the seirii-raajor axis in 
floating point inches. 

B is the length of the semi -minor axis in 
floating point inches. 

f\L'FHn is the angle of orientation (rotation) of 

the major axis in floating point degrees. 
This angle is measured from the horizontal 
to the major axis. Counter-clockwise 
Xs defined as the positive direction. 

THETAO is the starting angle between the radius 
vector (from center) and the major axis 
in floating point degrees. Counter- 
clockwise is positive. 



2.20.2 



TIIETAF is the final angle between radius vector 
and major axis in floating point degrees. 
Counter-clockwise is positive. 



IV 



For a closed curve TIIETAF = TIIETAO + 360.0 

=2 move to starting point with pen down 

(trace a line to starting point) . 
=3 move to starting point with pen up. 



Other Subroutines Required: 
Storage required: 
Restrictions: 

Original Source: 
In-house Submitter : 
In-hour.e Contact: 



PLOT 

DECsystem-lO (263) g Locations 

No ellipse is generated 

when , 

TIIETAO = THETAF 

Langley Research Center 

Tom Talbot/iaiL 

Tom Talbot/iaiL 



2.20.2.1 



r 

r A TF.«;T program to e;HOW THF \]^F OF «;tjRRO(IT INE FLIPS 

r 

''ALL PLOT«;(n,n,i6) 

r 

r CPT OPTGTN 
r 

(-All Of nT(^.,-n ,,_fl) 

CALL PLOT(.«>,?,5»-3) 

r 

C DRAW X AND Y AXF.«; 

r 

CALL AXIci?(0,,n,,3,,n,,i,,_4.,flHX - AXI c ,8,-l .XSC) 
CALL AXI<^?{n, ,o.,6.,0, ,1,,+6.,8HY - AXI c fe,-! ,YSC) 

r 

r noAw PLlTPcp 

r ■ 

CALL FLT P? (?,5,3. ♦!,,?,, 30. ,0,, 360,, 3) 

C 

C PART OF AN FLLTP.S CAN BF DRAWN BY CHANGING THF 6TH AND 7TH ARGUMFNT*; 

r 

r PND OF PLOT 

r 

CALL PL0T(ft,,ri,,-3) 

r 

r r;PMFPATF A rTPrilLAP FLl TP<^F 

r 

r CFT OPTGTM 

r 

CAf L OLOT(n,^_Ti,,_^) 

CALL PLOT(,s;,,'7,-3) 

r 

r DRAW X AND Y AXIS 
r 

CALL AXIS2(],.?.,3.,0,,1.,-4,,8HX - AXlS,8,-l ,XSC) 
CALL AXI«;2(1. ♦?,,5.,n.,l,,+6. ,8HY - AXI S«8,-l , YSC) 
r 

r DRAW FLLTPSF 
r 

CALL FLLTPS(3.,5.,l.,l,,45.,0,,36n,,3) 

r 

r FND PI OT 
r 

CALL PL0T(R.»n,,-3) 

STOP 
FMD 



2.20.2.2 



5-r- 



3-- 



»— • 

X 

CE 

I 



2 -- 








X - AXIS 



2.20.2.3 



5-T- 



4 -- 



3 — 



CO 

t— ) 

X 

CE 

I 

>- 



1 -- 




X - AXIS 



PICTORAL APPLICATIONS 
S306XX 



2.21 



Identification; Subroutine WORLD S30600 

Purpose 

To generate projections of the entire earth or 
selected portions of it as window plots. 

Usage ; FORTRAN 



where , 



Calling Sequence; 

CALL WORLD (KARD, lARGl , ARG2 , RDER, FILEN) 



KARD When KARD=1, a new plot is initiated 

according to the scaling factors input by 
the user. When KARD=2, a second plot is 
executed using the same scaling factors as 
as the previous call . 

lARGl Projection nuinber, specifies type of 

projection to be plotted. If no input 
number given, Mercator projections will be 
assigned by the program. 

ARG2 Set equal to 0. This argument is included 
here but intended only for possible future 
use . 

RDER Array containing 12 values 



(1) 
(2) 



(3) 



(4) 
(5) 
(6) 



(7) 



Central meridian ixxx.xx 

Floating point variable indicating 

spacing at which longitude lines are 

plotted XX. X 

Floating point variable indicating 

spacing at which latitude lines are 

plotted XX. X 

Set to 0. (No current use) 

Set to 0. {No current use) 

Floating point variable indicating 

northern latitude boundary of area to 

be plotted ixxx.xx 

Floating point variable indicating 

southern latitude boundary of area to 

be plotted ixxx.xx 



2.21.0 



(8) Floating point variable equal to eastern 
longitude boundary of area to be plotted 
±xxx.xx 

C9) Floating point variable equal to western 
longitude boundary of area to be plotted 

±XXX .XX 

(10) Set equal to (No current use) 

(11) Set equal to (No current use) 

(12) Integer variable l=grid, 0=no grid 

PILEN Data file indicator. 

FILEN(l) is data file WORLD. DAT for plotting 
the earth 

FILEN(2) is data file US. DAT used to draw 
boundaries of the continental United States 
(excluding Alaska) . 



Comments 



The map plot routine is written in FORTRM as a 
subroutine which will generate projections of the earth or 
selected portions . The user must initialize parameters , 
selecting the routines which will satisfy the user input 
requirements. The world plot and Continental U.S. boundary 
data is available on DECtape, 

There are four types of cylindrical projections 
available; 1) Mercator, 2) Equal Area, 3) Stereographic , and 
4) Equal Spaced. The axis of the cylinder is parallel 
to the axis of the generating globe. The basic features and 
any restrictions or limitations for each projection are 
presented below. The two digit number in parentheses which 
follows the projection name is the assigned nvmiber for that 
projection. 

2.21.0.1 



Mercator Projection (11) 
Basic Features 

1) Orthomorphic projection; i.e., shape preserved 
for small areas. 

2) Areas grossly exaggerated in polar regions. 

3) Windowing available. 

4) User may choose any longitude for central meridian. 
Restrictions and/or Limitations 

1) Latitude will not be plotted above 85 degrees for 
a full earth projection. 

21 Latitudes will not be plotted above 89 degrees 
north or 89 degrees south for a window projection. 

Equal ^rea (12) 
Basic Features 



1) Shape badly distorted (compressed) in polar 



regions. 



2) Equal -area property preserved over entire projection, 

3) Windowing available . 

4) User may specify any longitude to be central 
meridian. 

Restrictions and/or Limitations 
Latitude will not be plotted above 80 degrees 
north or below 80 degrees south for a full earth proJGCtion. 



2.2-1.0.2 



Cylindrical Stereographic (13) 
Basic Features 

1) Shape slightly deformed in polar regions. 

2) Scale most accurate between latitudes 45N and 
45S (latitudes where cylinder intersects globe) , 

3) Best compromise projection for shape, area, scale. 

4) Windowing available. 

5) User may specify any longitude to be central 
meridian. 

Equal Spaced. (14) 
Basic Features 

1) Equidistant projection-scale correct along all 
meridians and equator. 

2) Windowing available. 

3) User may select any longitude to be central meridian, 



Other Subroutine Required; 
Storage Required: 
Original Source: 
In-house Submitter: 
In-house Contact: 



PLOT, SYMBOL, NUMBER 



7162 



8 



S. Ross and Company 
Tom Talbot 
Tom Talbot 



2.21.0.3 



r 

r A tf<;t program to ^how thf ii<;f of worloplot 

r 

DIMFNSION POFPd?) ,PDFR2 ( 1? ) ,RDFR-? ( 12) ♦RDFR4 ( 12 ) .RDFRS ( 12 ) 
DTMFKKION FTLFN(?) 
INTFGFR FILPM 

DATA RDFR/-100,,?n,-,?n,,n.,0,,85.»-60,,180,,-180»»ft.»0,,l/ 
DATA ROFP2/-O«>,,?n,,?o,,n.,0.,5n,, + 2A,,-60,,-130,,A,,0,»0/ 
DATA RDFP'^/1S.,2n, ♦?n. ,0. ,0, ,"50. ♦?«>, ♦i»n, ,-1 0, ,0. ♦O. ,1/ 
DATA RDFPA/i^«>.,ir,,i'^.»n.,n. ♦5r).,'^n, ,155.,n'>,,0»,0.,l/ 
DATA RDFP5/l?5.5,,5.»0.,n, ♦2n,,-l«5, , 160. , 90. ,0, ,0» , 1 / 

r 

r vnRLD,DAT = IMP|fT FTLF FOP VORLD PLOTTING 
r U5^.DAT = INPMT FTLF FOP l\,<^, PLOTTING 

r 

DATA FILFN/'WORLD' ♦«(]<;♦/ 
TALL PL0T.<;(n,n,1ft) 

r 

r ^FT ORTGTN 

r 

CALL PLOTfO. ,-11 .»-'=^) 
CALL PL0T(n.,1 .♦-^) 

r 

r <^,FT PARAMFTFP*; FOP MFRCATOP PROJFCTION 
r 

APG?=0. 

IAPG1=11 

KAPn=i ^ 

r 

C PLOT rONTINFNTAL ll.«;. OVFRLAYFD WITH «;TATF BOUNDARIES 

r 

TALL VORLOfVARD^TARGI ,ARG2,RDFR?,FTLFN(1 ) ) 

(fAons*? 

TALL MOPLDdCARD, TAPGl , ARG? ,RnFR2 , PT LFN ( 2 ) ) 
r 

r FND PLOT 
r 

CALL DL0T(l«5,,n.,-?) 

(fARD=1 

r 

r WORLD PLOT - MPRCATOR PROJECTION 

r 

CALL W0RLn(lfAPn»TAPG1 ♦ ApG?,RnFR,FTLFN(l M 
GALL PLOT(l2.,n.,-^) 

r 

r WORLD PLOT - FOUAL ARFA PROJECTION 
r 

IAPG1=12 

CALL W0RLD<tfARD,IARG1 ,ARG2,RDPR,FTLFN(T )) 

CA| L PLOTd?.,/^,,--^) 
C 

C WORLD PLOT - CYLINDRICAL .<;TERFOGRAPHIC PROJECTION 
r 

IARG1 = T^ 

CALL WORLDfk'APD.TARGI ♦ARG2,RnFR»FILFNn) ) 



CALL PLOT(r?.,o.,-?) 

r 

r wn^LD PLOT - FnUALLY «^PACFD PROJECTION 

r 

IARG1=1A 

CALL WORLD(|CARr),IAPGl ♦ ARG2 ,RDFR,FILEN ( 1 )) 
CALL PLOT(l?.,n,,-^) 

r 

C MERCATOR PROJECTION PLOT«^ 0F9 
C MEDITERRANEAN AREA 

'" JADAM ANH. *'nPFfl 

r ea<;t tndtf^ 

r 

IARGl = n 

CALL WORLn(k'ARn, TARG1 ,ARG?,Rr)FR?,FTLFN(n ^ 
CALL DLOTn?.,n.,-^) 

CALL WOPLrXffARD, TAPG1 , ARCn? ♦RDFRA^FTl FN ( 1 ) ) 
CALL PLOT(l?,,n.,-^) 

CALL W0RLD{KARD,IARC1,ARG?,RDER5,IARG3,FILFNM)) 
CALL PL0T(1?,,0.,-?) 

«^TOP 
FMD 
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Identific ation; Subroutine SHA DE S 30601 
Pu rpose ; 

To shade a polygonal area with equidistant 
parallel lines. 

Usaqiii FORI RAN Calling Sequence 

CALL SHADE (XVERT , YVERT , NVERT , XWORK , YWORK, 

XTEMP, NTEMP, ANGLE, DIST) 



where , 



XVERT 

YVERT 
NVfRT 

xwIrk , 

YWOFK 
XTEMP 



V 



NTEMP 
ANGLE 

DIST 



An array containing the X-coordinates 
of the vertices of the polygon 
(ordered) . 

Same for the Y-Coordinaten. 

The number of vertices in the polygon 

Two working arrays dimensioned in the 
users program at (lWERT+1) . 

A temporary array used to hold the 
X-coordinates of the end-points of the 
line segments needed for each shade line, 
Should be dimensioned at tv;ice the 
maximum number of times a shade line can 
be expected to intersect the polygon. 

Dimension of XTEMP. 

The angle (with respect to the users 
X-AXIS) of the shade lines. 

The distance between. shade lines (user 
units) . 



The following comn^on area contains the- ,|)araraters used to 
convert from user coordinates to plotter coordinates. These 
values must be set by the user in his program in the same 



coniii'.on eye-". . 

i 
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Before plotting, the following transformation takes place: 
U= (X-XOFF) *XSF+UOFF 
V= (Y-YOFF) *ySF+VOFF 
COMMON/SKDCOM/XOFF , YOFF ,XSF , YSF , UOFF , VOFF 



Other Subroutines Required: 
Storage Required: 
Original Source: 
In-house Submitter : 
In-house Contact: 



PLOT, INTERX, ORDER, 

DEC system- 10 (620) g Locations 

Jim Steinberg/FIN 

Tom Talbot/KIIL 

Jim Steinberg/F.IN 
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r ^ TP«;T oponPAM to ^HOW thF MciF of «;nRPO|(TTNF <;HAnF 

r 

COMMON /«;HDCOM/X0FF,Y0FF,X.SF»Y<;F,tl0FF,V0FF 

DIMFNSION XWORK(IO) ,YW0RK(10) ,XTFMP ( ] nO ) ♦XARRAY( 10 ) , YARRAY ( 1 ) 

DATA XARPAY/l,,?,,4.,3.,1,,?.,n.,?,,l.,1.-/ 

DATA YAPPAY/1.,U,^.,5.,4.,3.,2.«>,?.,U5,1./ 

r 

CALL "PLOT<?('^, 0,16) 

r 

r f.c-j OPTCTN 

c 

CAl L PLOT^n. »-i i ,,-'^) 

CALL PLOT?.*;,. S,-^) 

r 

r DPAW APFA TO RF c;HADFD( POLYGONAL) 
C 

CALL PLOT(XARPAY(] ) ♦YARRAY^ n ,3) 

no 1 T=?,in 

CALL PLOT (XA PRAY (I ) ♦YARRAY ( I ) ,2) 
1 COMTTKMIF 
r 

C DRAW <;HAnF LTNF<; O.75 INCH APART, 45 DFCRFF<; TO HORIZONTAL 

r 

XCF=1, 
VCF=T, 

ItOFFsO, 
VOFF-0, ' 

xoFF=n. 

YOFFsO, 
C 

CALL «;HAnF(XARPAY,YAPRAY,9,XWORK,YWORK»XTFMP,100,45.,n.25) 

r 

r FND PLOT 

r 

CAl L PLOT(fl, ,n, ^--^^ ) 
e;TOD 

FMD 



IDENTIFICATION; Subroutine SHADIT S30602 

PURPOSE ; 

To draw bars and to shade these bars in different ways. 

USAGE: FORTRAN Calling Sequence 

Call SHADIT (Xl, Yl, X2 , Y2, DX, IPLT) 
where, 

XI Left handx coordinate (inches 
Yl Lower Y coordinate (inches) 
X2 Right hand X coordinate (inches) 
Y2 Upper Y coordinate (inches) 
DX Distance between shade lines 
IPLT Type of shading desired 

OTHER SUBROfjTINES REQUIRED ; Plot 

I 

STORAGE REQUIRED ; DECsystem-10 (1280) location 

COMMENTS ; 

Subroutine SHADIT immediately draws the rectangle defined 
by the points XI, Yl, X2, Y2. These coordinates are passed to 
the subroutine in CALCOMP compatible units (inches) . 
Subroutine SHADIT then proceeds to shade the rectangle, 
depending upon the value of variable IPLT. Not that when 
IPLT=5, no shading is done, i.e., only the bar is drawn 
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The following are examples of the five types of shading 
available. 





IPLT=1 IPLT=2 




• • • 

• • • 



IPLT=3 



IPLT=4 



IPLT=5 



Original Source: 
In-house Submitter; 
In-house Contact: 



Glenn Farrell, KHL 
Glenn Farrell, KHL 
Glenn Farrell, KHL 



r 

r A TFST PPOnRAM TO <;HOW THF USF OF .«;URR0UTINF .<;HADIT 

r 

CALL PLOT<;(n,n,ift) 

r 

r '-.PT ORTr,TM 

r 

CALL PLOT (.<=;♦, «;,-?) 
r 
r OPAW X AND V AyFc; 

r 

CALL AXI<;2(.5,,5,5,,n. ,i,,-6,,8HX - AX I <;»8 ♦ + ? ♦XSO 
CALL AXT<^?(.'5,,5,5.,n. ,1.,+4.,8HY - AXI <; , 8 ,-1 ,YSC) 

r 

r nPAV RAPc WTTH nyPFFPFNT CRO<:<; HATCHINCt 

r 

CALL CHAOTTf.?,.'^,! .^♦•=^.'5,.1 ,1 ) 

CALL «^HAnTT(1 .'i,,5,l.R,?.'i,.1 ♦?) 

CALl cHADTTf?. ♦.5,?,4,4.,.'^?,4) 

r 

r TMo Pt OT 
c 

CALL PLOT(T^. ♦n. ,-^V 

«iTOD 

FNO 
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4 -- 
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X -XflXIS 



4.00 
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UNIQUE PLOTS 
S307XX 
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Identification: Subroutine THREE 



S30700 



Purpose 

To plot a three dimensional perspective of the data 
of a two dimensional array, interpreting the value of the 
array as the 'z' component and the indices of the array as 
the 'x' and 'y' components. The program will also plot a 
function using subroutine EXTSCH with a user supplied 
function. 

Usage : FORTRAN 



Calling Sequence ; 

CALL THREE (A, AMAX, DISC, MAXDMX, MAXDMY, STARTX, 
SKIPX, DIMX, STARTY, SKIPY , DIMY , THETA, PHI, OPAQUE, 
BLOCK, ARRAY, SCALEX, SCALE Y, SCALEZ , Kl , K2 , K3) 



where, 



A Array to be plotted 

AMAX Maximum value of array or function to be 
plotted 

MAXDMX If plotting array, this is the maximum 
value. If not, set to 1 

MAXDMY If plotting array, this is the minimum 
value. If not, set to 1 

STARTX Starting X value 

SKIPX X increment, normally 1 

DIMX X dimension of array A 

STARTY Starting Y value 

SKIPY Y increment, normally 1 
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DIMY 

THETA 

PHI 



OPAQUE 

BLOCK 

ARRAY 



Y dimension of array A 

Angle from Z-axis to line of sight 

Angle in Z-Y plane between X-axis and line 
of sight 

If OPAQUE= . TRUE . hidden lines are not plotted. 
Not used at this time 

Logical variable. If .TRUE, the array 
A is to be plotted. If .FALSE, the user 
has supplied a function to be plotted in 
the form SUBROUTINE EXTSCH (DUMMY, IX, 
lY, Z) where DUMMY is not used 
IX = X value (integer) input 
lY = y value 
Z = resultant function value 



SCALEX Length in inches of plot X-axis 

SCALEY Length in inches of plot Y-axis 

SCALEZ Length in inches of plot Z-axis 

Kl Do not plot in Y direction if .TRUE. 

K2 Do not plot in X direction if .TRUE. 
K3 , Dummy . Not in current use . 
Other Subroutines Required: PLOT 
Storage Required: 21221g 

Original Source: Scott Moffatt/KHL 

In-house Submitters Tom Talbot/KHL 

In-house Contact: Tom Talbot/KHL 
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r 



r A tf<;t pponpAM to «-,nnw thf ik.f of «;imr?oi»TrMF thpff 

r 

LOr,TrAL y^ ."^^ ,1^^ , APRAY,R| OrK,OPAO|1F 

IMTFGrp cTAPTX,<;KTPX,r)TMy,«;TARTY,<^K'TPY,OrMV 
niMFMc TOM A( 1 ) 

TAIL t'I.OT«?(o,n,T A,) 

r 

r <;FT OPIATN 

r 

CAl L DLOT(n.,-T| .,-?) 
TAIL OLOT(n,,4,e;,_-^) 

<;ft dtmfn,ston<; of a arpay (omittfd in this FXAMPLF) 



r 
r 
r 



MAVDMYsl 

MAyOMY=1 

r 

r r,pT y AMf^ Y PAWrnFc 

r 

STARTX=1 
SKTPX=4 

niMx=ini 

«;taotv = 1 

«;ktpv=6 

DTMY=r1 m 

r 



r «;FT «;rALTNn RAMGF, maximum VALflF,ANn VTFWING AMr.LF 

r 

.«;rALFY=q, 

<;CALF7=7. 
AMAX = «in, 
PHI=f,0, 

THFTA=7'=i. 

r 

r MAifF PLOT OPAOtjFP nO NOT PLOT APPAYfPLOT FUNCTION) 

r 

OPAO|(r = ,TP!)F, 
APPAY=,FAl <;r, 

r 



r PLOT TN X niRFrTION IF lf? = .TPllF.,lN Y DIRFrTION IF K1 = .TRI)F, 
r 

K1=,TRUF. 
K?=,TRIIF. 

CALL THRFF(A»AMAX,DUMMY,MAXDMX»MAXOMY»STARTX.SKlPX»DlMX*STARTY,.SKT 

*PY ^D IMY, THFTA, PH N OP AOl)F,RLOCK,ARR A Y,5;CALFX,,<;CALFY,5;CALF7 ♦<!♦!<:?, 

*ni(MMY) 

r 

r Fwn PLOT 

r " 

CALL PLOTnft.,n.,_^) 

<;top 

FNO 
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SPECIALTY ROUTINES 
S308XX 
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Identification ! Program DECKS S30800 

Purpose ; 

To plot program deck structures according to user specification 
primarily for use in documentation. The subroutine will plot 
cards and/or decks of variable length with printing space for 
40 characters. 

Usage; FORTRAN 

Calling Sequence 

Call Decks 

Comments ; 

The contents and spacing of the plotted decks or cards aire 
specified by the use of control cards in the user's data 
file. The control cards are identified by an asterisk 
in column 1. All non-control cards are copied on the 
plotter. An explanation of the control cards follows; 

*CARD Used to add an additional offset space for 
the next card to be drawn, or to end the 
drawing of a deck, 

♦COLUMNS To draw column numbers on next plotted card 
from the column position specified by the 
integer ending in colvimn 12 to the position 
specified by the integer ending in column 16 
More than one *COLUMNS card may be used to 
specify the printing of any column or 
columns . 

*DECK Indicates that a deck representation is to 
be drawn. The first card after the *DECK 
card that is not a *CARD or *COLUMNS card 
will be used as the front of the drawn deck. 
A blank card my be used if the front of the 
drawn deck is to be left blank. 
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All cards after the first card v/ill be used 
to draw labels on the top of the deck. The 
size of the drciv/n deck is controlled by the 
number of blank cards . 

A drawn deck is ended by a *CARD. 

Begin a new plot . 

Used to end the drawing of a deck setup. 
Other Subroutines Required: PLOT, SYMBOL, NUMBER 
Storage Required: 1441, 



*NEXT 
*STOP 



Restrictions: 

Original Source: 
In-house Submitter: 
In-house Contact: 



'8 

The user must take care not to 
draw a deck longer than the 
plotter paper size 

Langley Research Center 

Tom Talbot/KHL 

Tom Talbot/KHL 
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r 

r A tf«;t progpam to ."^how thf osf of subroutine decks 

r 

rALL PLnT«;(n,n,ift) 

r 

rA| L nFrv<; 

r 

r FNH PLOT 

r 

TALL t>l.OT< 1?. ,0, ,-?) 

«;Tnp 

FMD 

r 

r SAMPLF DATA FTLF^ TO RF IJSFD WITH SURROOTINF PFCKS 

r 

r FTLF MAMr HFK'.nAT 

r 

r THIS OFriC <^TRurTIIPF WAS ItSFO TO PROVIDF THF FXAMPLF 

r 

•RSFOIfFNTP 

*JOP 

♦PASSWORD 

♦rAPo 

<; FORTRAN 

♦ DFr< 
FORTRAN PPOnPAM 

♦TAPn 
sr>ATA 

♦COLUMNS 1 5 

#rO| MMNC TO TO 

♦ rni tiMN<; i^ is 

♦ HPrtf 

71. ft Q in 

nFri<? 

ACAPO 

♦ CARD 

FNin-OF-FILP 
♦CARD 

♦ CAPO 

♦ STOP 
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IDENTIFICATION: Proqram DOCPRT S30801 

Purpose; 

This program is used to reformat computer listings to 
page size. The file name and page numbers are placed on 
each page making it an ideal method for documenting a 
computer program or data files. 

Usage ; FORTRAN 
Comments: 



The procedure for executing this program is as follows; 

• MOUNT DTA: SAV/VID: 611/WLJ 

• R PIPj 

* DSK:/X = DTA_: DOCPRT. SAV,^ 

* +C|^ 



1. .RUN DOCPRT 

2. Type in a page number prefix. This is a one (1) 
character letter associated with the appendix into 
which these listings are to be placed. To reduce 
typing, xeroxing, and collating, all along listings 
in bound reports will be put in appendices numbered 
A-1, A-2, etc., for Appendix A; and B-1, B-2, etc., 
for Appendix B, etc. If the work "NONE" is typed, 
the listings will be produced without page numbers. 
This may be desired for some short listings to be 
inserted in the text of a report or in an external 
memo. This program can also be used for listing 
input files or control files, and appendix page num- 
bers may not be desired for these. 
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3. Type in the starting page number if needed. 

4. Type in a project or program name. This is 
printed on the top of each page of the listing 
surrounded by asterisks. 

5. Type in the full filename. A look up error 
occurs when the named file is not found. When 
this happens the program simply loops back to 
accept another filename. Type "DONE" to stop. 

6. When the program asks for "NEW APPENDIX PREFIX" 

or "NEW TITLE? P/T/<CR>:", a P returns the program 
to step 2 above, a T to step 4, and a <CR> to step 
5. Previously used titles and appendix prefixes 
with continuing page numbers are used if new ones 
are not asked for and accepted. 

7. After the program exists following a DONE file name, 
the documentation listings are in a file named 
DOCUI^.DAT. Print with /COPIES: as many copies as 
are i|ecessary to eliminate any xeroxing of listings. 

8. DISMOUNT DTA_:/Re4 

Other Subroutines Required; None 
Storage Required: 

Restrictions; 

During the execution your file is unamed to PTEMP.DAT. 
Therefore if for some reason the program does not do a 
normal exit look for your last file to be named PTEMP.DAT. 

Original Source ; Dr. R. Wassmutm, T. Talbot/KHL 
In-House Contact: Dr. R. Wassmuth/KHL 
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M OUNT riTA:SAVE/UID:611/tJ L 
Reouest Gueued 
Uaitinsi* * *2 '"C's to Exit 

ZSL 
.'in 

♦ SYS ^S- 

27 :i001rA7 7 TTY4;? SYSTAT 12fSPY RN * 

* means tiKecute Only 

RES^ 
pfT2 » PC02 r PT32 » PT 1 y DSK A » DSKB r DSKD f PVT6 » DSKC r TRAS f PVT3 f PCOl » BIZ4 f PVT7 » UPAS 
DPA7 f DPB6 f CriROOO » LPTOOO » 00 1 » 002 ? DTAOOl » 002 r 003 » 005 r MTAOOl r 004 » PLTOOO 
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f^77 TTY43 SYSTAT :l2fSPY RN * 

$ means Execute Onlvj 



Busy dev.ices{ 




Device Jab Why 


Logical 


nTA004 27 3S 


SAVE 


♦R PIPJ 




*DSKB : /X:^=DTA004 : DOCPRT . 


► SAM^ 


*2£ 




♦RUN DOCPRT. SAV 2 





%runn:[ng ka-10 code on a ki-io 

input appendix prefix— none for no page numbers? q,^ 

type in starting page number; ^^ 

input project or program title? 
example hydras system programs 

supplementary plotting routines ^ 

input full filename—done to stop? tencur.fip ^ 

lexceeding quota on dskb3 

new appendix prefix or new title? p/t/<cr>:il 
input full filename— done to stop: done ^ 



CPU time: 10.93 ELAPSED TIME? 2:33.85 
NO EXECUTION ERRORS DETECTED 



.DIR IiOCUM.DAT 

l.OCUM DAT 287 «)57> 31-hBi 77 USKEC 1 3001 » 477:1 



.PRINT DOCUM.DAT 

Total of 287 blocks in 1 file in LPT reouest 

♦PRINT 

The oueue is enriptw 

♦ DEL riOCUMlDAT 

Files deleted* 

COCUM.DAT 

287 BlocM.s freed 
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MISCELLANEOUS 
S309XX 
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SECTION 3 



3.0 



PLOTTING PROCEDURES 

the riecsystem-lo has been checked out with the new 
calcomp 915 stand alone controller/plotter for both the 7 or 9 
track magnetic tape drives* 

the following procedure should be used when pen plotting 
or film plotting on magnetic tapes generated on the decsystem-10* 
<cr> means carriage return, 

!♦ log onto the decsystem~10 computer* 

2. use the resource command to check mta drives avail- 
able where mta0»mtal»mta2»mra3 are 7 track and mta4f mta5»mta6f 
mta7 are 9 track drives. 

example: 

.resources<cr> 

3. request (by use of the mount command) a scratch magnetic 
plot tape on a mta unit or if desired* your own plot tape 
number. either a 7 or 9 track drive may be used depending on 
what drives are available from the resource command. 

7 track is the system default if track is not specified in 
the mount command. 

examples: 

.MOUNT MTA: 16/REELID:1234/VID:'J0NES-7 TRACK PLOT TAPEVWENABL<CR> 
.MOUNT MTAJ 16/REELIDt 1234/VIDJ ' J0NES~9 TRACK PLOT TAPEVWEN<CR> 

WHERE 16 IS A LOGICAL UNIT NUMBER r 1234 IS THE MAGNETIC 

TAPE NUMBER* AND JONES IS YOUR NAME. A SCRATCH TAPE 

IS PREFERRED BY THE OPERATOR BUT YOUR OWN TAPE MAY BE USED IF DESIRED. 

IMPORTANT~BE SURE TO PUT APOSTROPHES AROUND YOUR MESSAGE ON VID SWITCH. 

4. SET BLOCK SIZE 82 FOR FORTRAN-4 FOROTS ONLY* NOT FORTRAN-10 
EXAMPLE? 

.SET BLOCKSIZE 16: 82<CR> FOR FOROTS F~4 ONLY 

5A. FOR PEN PLOTS ~ EXECUTE YOUR PLOT PROGRAM USING THE STANDARD CALLS 

available in the calcomp manual for pen plotters sept. 1969 for the 
ibm 7094. these subroutines include axis* line* number* scale* symbol* 
plot* buff* factor* newpen* where* and plots. 

example: 

.set blocksize 16 82<cr> for f-4 only 
.exec/f40/l0ader/f0r0ts pgm.f4*sys:94plt/lib<cr> 

where pgm is your program for plotting, 
till 3/31/77 when fortran iv will be off the sys area 
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FOR FORTRAN-IO: 

.EX PGM.F10»NEU}94PLT1.REL/SEA<CR> 

SB. FOR PEN DRAFTING OR LOWER CASE LETTERING |^J"^f §^J^. J^f 
PLOT PROGRAM USING THE STANDARD CALLS IN THE CALCMPGR^^^ 

[^??^^.?^f sirr/u^^ T.T\.:z.T'v.Ts sfL fx^^L^A^Nf T^r 

^%Y.'.l%uZ.ir.%^. ;^rs?oSyTc^A^^fs?HJrVf NEU VERSION., ^ 
PLS?BiFF FACTOR, NEWPEN, WHERE, PLOTS, SYMS <NEU)»SCURV <NEW)'1-LEX (NEW), 
AND LBYTE (NEW ASSEMBLY LANGUAGE ROUTINE RbQUIRED BY SYMBOL). 

example: 

.SET BLOCKS I ZE 16 82<CR> FOR FORTRAN --4 ONLY 
. EXEC/F40/L0ADER/F0R0TS PGM . F4 , SYS J DRAFT/LI B<CR> 
TILL 3/31/77 WHEN FORTRAN IV WILL BE REMOVED FROM THE SYS AREA. 
FOR F0RTRAN~10: 

.EX pgm.fio,new:drafti.rel/sea<cr> 

5C. FOR FILM PLOTS ~ EXECUTE YOUR PROGRAM USING THE STANDARD 

CALLS AVAILABLE FOR FILM PLOTTERS (MAR. J?!? MANUAL) FOR THE !»« f ^94 . 

THESE SUBROUTINES INCLUDE - AXIS, LINE, SCALE, NUMBER, SYMBOL, LALCMh , 

reset, bufout, and buff, 
example: 

.SET BLOCKSIZE 16 82<CR:> FOR FORTRAN-4 ONLY 
. EXEC/F4O/L0ADER/F0R0TS PGM . F4 , S YS : 835PT<CR> 
TILL 3/31/77 WHEN FORTRAN IV WILL BE REMOVED FROM SYS: 
FOR FORTRAN-IO: 
.EX PGM.F10,NEW:835PT1.REL/SEA<CR> 

6. DISMOUNT YOUR MAGNETIC TAPE WHEN PLOTTING JOB IS COMPLETED. 

EXAMPLE : 

.DISMOUNT 16:<CR> 

7. CALL THE I/O DESK AT X42674 TO PLOT YOUR TAPE. 

FOR F='EN plots: 

EXAMPLE OF WHAT IS NEEDED: 

PLEASE PLOT 7 TRACK TAPE NO. 1234, JONES, ACCT. E61616A, 
12- PAPER, BALLPOINT PEN, 5 PLOTS, X41121, TAPE TO BIN E56. 
YOU MAY ALSO FILL OUT THE PLOT REQUEST AT I/O ROOM 140-C. 

FOR FILM plots: 

EXAMPLE OF WHAT IS NEEDED: __^ ^^^^ ^^^z.^a 

please PLOT 7 TRACK TAPE NO. 1234, JONES, ACCT. E61616A, 

35 MM FILM, 5 PLOTS, X41121, TAPE TO LIBRARY. 

YOU MAY ALSO FILL OUT THE PLOT REQUEST AT I/O ROOM 140-C. 

3.2.1 



important: your call to plots or CALCMP should have a 16 IN THE 
3RD ARGUMENT. 

EXAMPLE— CALL PL0TS<0»0f 1A»0) FOR THE PEN PLOTTER 

CALL CALCMP(0r0rl6,0) FOR THE FILM PLOTTER 

PREVIOUS DECSYSTEM-10 PLOT PROGRAMS WILL HAVE TO BE CORRECTED AND 
RECOMPILED TO MATCH THE CALLS ON THE 915 CONTROLLER. 

FOR THOSE USING KHLIB OR OTHER SPECIAL PLOT PACKAGESf LOAD WITH 
LINK BEFORE EXECUTION TO AVOID GLOBAL ERRORS. 

EXAMPLE 

.R LINK<CR> 
*MYPLT<CR> 
J!?JSYSJKHLIB/SEA<CR> 
#SYSJDRAFT/SEA<CR> 
*/NOSYSLIB<CR> 

*/GO<CR> 

NOTE THAT THIS WILL LOAD IN ONE COPY OF ALL MODULES NEEDED AND 

THAT /SEA (SEARCH SWITCH) IS USED FOR LINK LOADING. 

THE /LIB SWITCH IS USED FOR LOADING WITH /LOADER. 

IF GLOBAL ERROR OCCURS r RESTATE THE LIBRARY TWICE WITH 

THE MODULE THAT IS MISSING. 



3.2.2 



TWO PLOTTING CAPABILITIES EXIST ON THE DECSYSTEM-10 FOR PLUI.xixu 
ON THE TEKTRONIX SCREENS: 

1. CALCOMP PREVIEW PACKAGE WHICH ALLOWS A USER TO PREVIEW 
A CALCOMP PLOT ON THE TEKTRONIX SCREEN. 

•?♦ ADVANCED GRAPHICS PACKAGE WHICH ALLOWS A USER TO WRITE 
HIS OWN TEKTRONIXS GRAPH PACKAGE USING THE TEKTRONIX 
SUBROUTINES. 



procedure: 



1. LOCATE A TEKTRONIXS TERMINAL AVAILABLE FOR USE. CALL 
RUSS NAHIGIAN AT X42015 OR SEE HIM IN BUILDING 3 ROOM 127 
FOR LOCATION OF TERMINALS AND AVAILABILITY OF USAGE 

AND FOR PLOTTING HELP. ALSO, TEKTRONIX MANUALS AND 
DOCUMENTATION ARE AVAILABLE. 

2. TURN ON THE MODEM ATTACHED TO THE TETRONIX SCREEN. 

3. TURN ON THE TEKTRONIXS UNIT (SWITCH IS UNDER THE KEYBOARD 
ON THE RIGHT SIDE OF THE FRONT PANEL). 

4. NOTE THE BAUD RATE OF THE TEKTRONIX UNIT. ^T MAY BE 300 
BAUD OR 1200 BAUD. THE BAUD RATE SWITCH IS ON THE BACK Oh 
THE TEKTRONIXS UNIT NEAR THE OUTGOING CABLES. DO NOT CHANGE 
THE SETTING BUT DIAL IN THE PROPER BAUD RATE PHONb NUMBER. 



D. 



DIAL THE APPROPRIATE PHONE NUMBER FOR THE BAUD RATE AND 
ATTACH THE PHONE TO THE MODEM. (IN SOME CASES A WHITE 
BUTTON ON THE PHONE MAY HAVE TO BE ACTIVATED BY 
LIFTING). THE MODUM LIGHT AND TEKTRONIX UNIT SHOULD 
BOTH INDICATE AN ON-LINE LIGHT. 

6. CLEAR THE TEKTRONIX SCREEN WITH CLEAR SCREEN 
CHARACTER ON THE KEYBOARD. (AFTER WARMUP PERIOD). 

7. TURN ON THE HARDCOPY UNIT UNIT IF DESIRED (AND AVAILABLE). 

8. CTRL C THE TERMINAL TO INSURE CORRECT CONNECTION, CLEAR THE 
INPUT BUFFER, AND RETURN THE PERIOD. 

9. LOG ONTO THE DEC -10 COMPUTER. 

10. INITIALIZE THE TELETYPE FOR CALCOMP PREVIEW. 

A.) .ASSIGN TTY 16 <CR> 

B.) .TTY FORM <CR> 

C.) .TTY NO CRLF <CR> 

D.) .TTY NO FILL <CR> 

WHERE <CR> MEANS THE CARRIAGE RETURN KEY 

11. EXECUTE YOUR PLOT PROGRAM. 
A.) FOR CALCOMP PREVIEW 

example: .EXEC/F40/LOADER/FOROTS MYPLOTf SYSITEKPLT.REL/LIB <CR> 
OR 
TILL 3/31/77 WHEN FORTRAN IV WILL BE REMOVED FROM SYS: 



\\ 



B.) FOR ADVANCED GRAPHICS USAGE 
example: ♦EXEC/F40/L0ADER/F0R0TS MYPL0T»SYS:TEKADG.REL/LIB <CR> 
TILL 3/31/77 UHEN FORTRAN IV WILL BE REMOVED FROM SYS: 

WHERE MYPLOT IS YOUR PLOTTING SOFTWARE PACKAGE. 

12. THE SCREEN WILL BLIP AND ASK YOU FOR OPTIONS? RESPONSES ARE: 

A.) 1 <CR> PLOT MY COMPLETE GRAPH. 

B.) 2 <CR> THE USER IS ALLOWED TO RESET THE SCREEN WITH 
NEGATIVE PEN OPTION. 

C.) 3 <CR> RESET ORIGIN OF EACH PLOT AND RESET SCREEN WITH 
NEGATIVE PEN OPTION. 

D.) 4 <CR> NOT USEDf BUT ALLOWS THE USER TO SUPPLY HIS OWN 
OPTIONS. 

13. THE TERMINAL WILL AWAIT ANOTHER RESPONSE. 
A.) C CONTINUEfPROCEED WITH PLOT 

B.) E ERASE THE SCREEN. 

C.) H HARDCOPY OF SCREEN IMAGE. 

D.) Q STOP PROGRAM 

E.) S SKIP NEGATIVE PEN CALLS. 

A QUESTION FOLLOWS ASKING HOW MANY SKIPS? 
example: 1 <CR> ONE SKIP 

E.) W WINDOW ADJUST. ALLOWS USER TO RESET THE PLOT 
ORIGIN. TWO QUESTIONS FOLLOW. 

1.) COORDINATES OF THE ORIGIN? (INCHES) 

example: -1»-1 <CR> (best start origin POINT) 

2.) WIDTH AND HEIGHT OF THE PLOT? (INCHES) 
example: A»8.5 <CR> 

3.) TYPE C <CR> TO CONTINUE AFTER COORD. ARE ENTERED. 

F.) R RESET WHICH ALLOWS THE USER TO SELECT ORIGINS 

OTHER THAN THE LOWER LEFT HAND CORNER OF THE SCREEN. 

QUESTIONS FOLLOW ASKING X AND Y COORDINATES OF THE 

ORIGIN? 

example: 2»3 <CR> 

? THIS MESSAGE 

14. YOUR PROGRAM SHOULD PLOT ON THE TEKTRONIX SCREEN. 

15. MAKE A HARD COPY OF THE PLOT IF DESIRED. 

16. PLOT MORE GRAPHS (AFTER RESETTING THE SCREEN) OR LOG OUT. 



-* «\ 



ONE ADDITIONAL NOTE J 

SOME USERS WANT TO USE SPECIAL PLOT PACKAGES qiirw ao kui rx. 
LOAD THESE SPECIAL ROUTINES WITH LINK TO aCoiS SlObL ERROrL'""''^ 

EXAMPLE 



R LINK<CR> 
*MYPLT<CR> 

*Sys:KHLIB/SEA<CR> 
*SYSJTEKPLT<CR> 
*/N0SYSLIB<CR> 
«/G0 



THE /SEA (SEARCH SWITCH) IS ONLY USED WITH I TWk- 
THE /LIB SWITCH IS USED FOR LOADING Si ¥h /LOADER. 

A^tSJs tLe^^'^" ""^ '""^ '"^"^ '='''^'^'^«^" ^S ^0^ A^^AILABLE 



3.2.5 



